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CHAPTER 1 


INTRODUCTION N72-27674 

J.C. Hung 


This is a technical report on a study of strapdown platform technology. 
The study was initiated by the National Aeronautics and Space Administration, 
George C. Marshall Space Flight Center, and was performed by the University 
of Tennessee under Contract NAS8-27296 over a period of one year, beginning 
April 16, 1971. 

Strapdown Principle 

Strapdown navigation platforms have emerged as a new generation of 
navigation equipment. The basic principle involves 1) the mounting of sev- 
eral gyros and accelerometers directly on the vehicle body and 2) performing 
the equivalent function of a stabilized platform by maintaining the coordi- 
nate reference in a computer.^ Three-axis body angular rates can be sensed 
by a set of gyros and/or accelerometers. This information can then be inte- 
grated over a short period of time to provide angular increments in the body's 
coordinate system. If the increments are sufficiently small and if the as- 
sociated digital computer can perform sufficiently rapid computations, then 

2 

an accurate angular reference can be maintained. 

There are several important advantages in using a strapdown system, 
as compared to a stabilized system. First, the nonexistence of a gimbaled 
platform eliminates all errors associated with the platform stabilization. 
Secondly, being free from gimbal lock allows an all attitude operation. 
Thirdly, because of the absence of mechanical platform gimbals, the system 
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can be smaller in size, lighter, more rugged, and easier to maintain. The 
last advantage, but probably the most important one, is the ease of the 
strapdown system to provide redundancy, which allows for an increase in 
system reliability. These advantages are achieved at the expense of a 
more complex electronic unit and a more stringent computer requirement. 
However, the progress of integrated circuit technology has made this trade- 
off worthwhile. 

In general, the development of strapdown technology can be divided 
into two levels, namely, the platform system level and the sensor subsys- 
tem level. At the platform system level, the main consideration is how to 
use a set of sensors in a manner which will yield the maximum hardware re- 
liability and give the best information accuracy. At the sensor subsystem 
level, the ultimate goal is to achieve the most accurate and reliable sensor, 
which depends on the sensor’s principle of operation, its mechanical con- 
struction, and the associated electronic circuitry. 

Rebalance Loop 

The crux of a strapdown sensor system is the use of rebalance torques 
to continuously restore the position of sensor output axes. The rebalance 
feature results in a much simpler computation requirement, reduces the ef- 
fect of nonlinearity, and avoids the occurrence of gimbal lock. Fig. 1-1 
is a schematic of a rebalance gyro loop. An input body rate causes a pre- 
cession about the gyro output axis. Motion of the output axis is sensed by 
a pickoff which passes this information to an electronic unit. A train of 
current pulses is generated by the electronic unit and the shape of these 
pulses is determined by the input signal. The sign of each pulse is deter- 
mined by the sign of the angle of gyro output axis while the amplitude and 
the width of each pulse is accurately maintained. Thus, each pulse repre- 
sents an incremental change of the output angle. These pulses are used for 
two different but related purposes: 1) as the current input to the gyro 

torquer to rebalance the output axis, and 2) as the discrete-data input to 
the computer for maintaining an attitude reference in the computer memory. 




4 


There are four major requirements which dictate the design of a re- 
balance loop. The first is the stiffness and the stability of the rebalance 
loop, which governs the dynamic error introduced by the sensor. The loop 
must provide a sufficiently large maximum torque capable of rebalancing the 
largest probable sustained angular rate or linear acceleration input. Thus, 
it calls for a well compensated high gain feedback loop. The second require- 
ment is placed on the characteristics of the generated current pulses. To 
achieve a good output resolution, a small pulse width ii necessary. However, 
the strength of each pulse must be sufficiently large to produce an effective 
rebalance torque. Third, the pulse frequency must be capable with the fre- 
quency response characteristics of the torquer to obtain good linearity over 
its entire dynamic range. The fourth requirement is that the shape of each 
current pulse must be sufficiently insensitive to the unavoidable variations 
of electronic components to avoid the "scale factor error". This assumes 
that the signal which the computer receives from the sensor will match the 
rebalanced mechanical angular motion of the sensor output axis. By an intel- 
ligent choice of rebalance technique, a careful design of electronic circuitry, 
and a properly compensated rebalance loop, this requirement can be achieved. 

High Reliability via Redundancy 

One of the significant advantages of a strapdown system is the ease 
with which sensor redundancy may be implemented. The system must be capable 
of sensor performance management which includes the functions of failure de- 
tection, faulty sensor identification, system reconfiguration and, possibly, 

sensor recalibration. A number of sensor redundancy configurations h«-.e been 

3—6 

developed using simple— degree-of— freedom sensors. Continued research in 

this area is still needed to improve the system reliability and accuracy. In 
principle, two-degree-of-freedom gyros can also be used in a strapdown system 
where each gyro will provide attitude information along two coordinate axes. 

It is important to know to implement this class of strapdown systems and how 
they will compare with other classes of systems in performance, reliability, 
and complexity. 
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Attitude Algorithms 

Two main difficulties of a strapdown navigation system, not shared 
by a stabilized navigation system, are the scale factor errors and the co- 
ordinate transformation errors. The former was mentioned in the previous 
section. Three factors which cause coordinate transformation errors are: 

1) insufficient speed of computation, 2) truncation errors, and 3) round- 
off errors. The most obvious approach would be to use a fast and large 
computer; but these errors can also be reduced by the use of a more effi- 
cient computation algorithm. Several transformation algorithms are known; ^ ^ 
however, at the present there are no guidelines which can be used for the 
selection of the best algorithm. 

Rebalance Electronics 

To have an accurate and reliable overall system, it is necessary that 
the rebalancing torquer be driven by very precise current pulses of well de- 
fined duration, amplitude, and shape. Therefore, the electronic pulse gen- 
erator must provide such an output over the range of environment expected. 
Because of the stringent performance demanded of the electronic pulse gener- 
ator, well designed electronic circuitry, which makes use of the best elec- 
tronic components and novel circuit concepts, is needed. 

Report Organization 

Chapter 2 of this report presents the analysis of the pulse rebalance 
loop. Chapter 3 is devoted to the development of a new concept of high re- 
liability strapdown systems using two-degree-of-freedom gyros for attitude 
sensing. Chapter 4 reviews the foundation of various transformation algo- 
rithms and presents their realizations by digital differential analyzers. 
Chapter 5 contains results of rebalance electronics analysis. Chapter 6, the 
last chapter, proposes problems for further study. A reference section is 
included at the end of each chapter except Chapters 1 and 6. 

Related Publications 

Several other publications reporting the partial results of the study 
conducted under this contract are listed below. 
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CHAPTER 2 


N72-27675 

REBALANCE LOOP ANALYSIS * 


J.C. Hung and W.L. Rowe 


ABSTRACT 


This chapter presents the analysis of the 
pulse rebalance loop where the pulses are of con- 
stant amplitude type. Properties of such systems 
will be discussed. Sensor errors caused by the 
distortion of the pulse shape, the dc bias, and 
the sampling process will be analyzed. 


INTRODUCTION 

The function of the sensor rebalance loop is to continuously restore 
the sensor axis from any variations to its aligned position. These variations 
are caused by changing attitude conditions of the vehicle. By properly mea- 
suring the rebalancing signal, the vehicle attitude can be determined. 

While several different types of signals may be used to rebalance the 
sensor, the needed measurement accuracy of the attitude change requires a 
calibrated digital signal in such a form that each pulse has a definite weight 
associated with it. Such a rebalance loop is shown in Fig. 2-1. The purpose 
of the pulse generator is to analyze the continuous signal from the pickoff 
and provide a number of pulses to rebalance the sensor and at the same time 
provide highly calibrated information for the computer. 

According to the form of rebalancing pulse stream, there are several 
different pulse-rebalance techniques as shown in Fig. 2-2. The binary tech- 
nique has pulses of constant frequency and amplitude. Essentially, the problem 
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posed to the pulse generator is the choice of proper polarity of the pulse. 
The width-modulated binary technique follows the same concept, yet adds to 
the pulse generator the problem of determining the width of the pulse with- 
in certain limits. The variable pulse width also creates problems in pro- 
viding the computer with quantized information since a simple counting tech- 
nique would not differentiate between pulses of like polarity but different 
widths. The ternary technique is casual, with a pulse occurring only when 
the pickoff amplitude exceeds some value. Both constant and variable width 
pulses can be used in the ternary method. Obviously, as the function of the 
pulse generator becomes more sophisticated, the circuit itself becomes larger 
and more complex. 

The purpose of this chapter is the analysis of possible errors in 
the strapdown sensor loop. Various errors due to the pulse generator are 
examined. Since the pulse generator is not ideal, one must consider errors 
such as the bias signal generated by the electronics and errors caused by 
the nonzero rise time of the torquing pulses. Knowing these errors and their 
effects, possible schemes of improvement are then considered. The effects 
of bias and of the shape of current pulses will be analyzed separately in the 
following sections. 


EFFECT OF BIAS ERROR 

Bias Error 

To produce the torque needed to rebalance the sensor, an electronic 
network is used to generate current pulses to drive the torquer, which is 
essentially an electromagnetic device. Ideally, the properties of this net- 
work include (1) a linear input-output relationship and (2) a precise cali- 
brate n. In order to have these properties, it is necessary that the shape 
of the current pulses be maintained and the net output current be zero when 
the input is zero. 

Effect of Bias Error on Sensor Accuracy 


From a knowledge of the system, which contains an integrator, it 
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Is apparent that the error due to bias drift Is cumulative and will re- 
quire a number of pulses per unit time to null output variations. The 
accumulation of these error pulses occurs in the computer which records 
all pulses as attitude variations seen by the sensor. It is the purpose 
of this section to analyze this effect and to define the period of these 
error pulses, t Q , as a function of the drift magnitude. 

The general system for the analysis of bias effect can be repre- 
sented by Fig. 2-3, where the bias is assumed to have a constant magnitude. 

A. Two types of analysis will be discussed in the solution for t Q . If 
the electronic network or pickoff has a threshold, T , a pulse will occur 
whenever the response of the sensor due to the bias drift becomes equal 
to the threshold. A second method of analysis, is to assume zero input 
condition with bias drift, such that the effective area under the output 

axis waveform is zero at t . 

o 

Method One . Given an initial response equal to the threshold, a 
pulse will occur which will torque the magnitude of the sensor response 
back within the limits of the detector, as shown in Fig. 2-4. Therefore, 
the time when the response due to the pulse and the bias drift again equals 
T g corresponds to t Q . That is, the time t Q satisfying the following equation 
is the period of the error pulses due to bias error. 

T + A * g (t) + B [u(t) — u(t — t ) ] * g (t) = T fi , (2-1) 

* O O 

where "*" represents the real convolution operation and B is the pulse height. 
Eq. (2-1) may be reduced to 

A * g (t) - B [u(t) - u(t - t )] * g (t) * 0 , (2-2) 

g g 


The equation describing the threshold case can now be reduced to 


A b (1 - e~ 8T ) 

2 " % \ 
as s (s + a) 



» 


0 


(2-3) 
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Figure 2-4. Quiescent Sensor Output. 
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which in the time domain yields 



B . B -a(t- 
a a 


"] u 


(t - T) - 0 


t**t 


(2-4) 


Neglecting the transient terms and writing Eq. (2-5) at time t Q , 



B 

a 


t 

o 



t) 


0 


(2-5) 


Solving Eq. (2-6) for t yields 


t 

o 


B T 
A 


(2-6) 


Assuming that the torquing pulse does not drive the sensor past the 
opposite threshold of the detector, the period of the error pulse, t Q , is 
given by (2-6). Even if the system is torqued past the opposite threshold, 
a net summation of positive and negative pulses will yield one net error 
pulse every t seconds. An illustration of the relation of bias amplitude, 

A, the period of the error pulses, t Q , and the frequency of the error pulses, 
N q , is shown in Table 2-1 and Fig. 2-5 for a pulse width and system time con- 
stant of one millisecond. 

Method Two . The second method of analysis involves observing the 
result of the bias drift and the pulse inputs to the system. If the pulse 
is applied at a time such that the sum of the positive and negative areas 
is zero, then the time when the cancellation occurs corresponds to t Q . Since 
the output at t is zero, the net effects of the bias drift and pulse inputs 
at this time are also zero. Fig. 2-6 illustrates this idea. 
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TABLE 2-1 

EFFECT OF BIAS ERROR ON PULSE OUTPUT OF THE SYSTEM 


Drift Magnitude 

Period of Error 

Frequency of 

Pulse Height, 

Pulses, t , 

Error Pulses, 

A/ J B| , units 

seconds 

N q , pulses/second 

10' 6 

10 3 

io ” 3 

2 x 10” 6 

500 

2 x 10~ 3 

A x 10” 6 

250 

A x 10~ 3 

8 x 10~ 6 

125 

8 x 10" 3 

io " 5 

100 

io - 2 

2 x 10" 5 

50 

2 x 10~ 2 


Note: a = 10 3 , t - 10 3 




A , Bias Drift, Amperes 


Figure 2-5. Period of Error Pulses Due to Bias Drift 


Sensor 

Output 



Figure 2-6. Equal Area Sensor Response. 
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The steady state response generated by the bias drift, y^(t) can be 


expressed as 


y l (t) " [ (s + a) | s-0 ] C » 


(2-7) 


or at time t 


At 

yi<t> - 


(2-8) 


Also, the steady state response due to the torquing pulse, y^(t), is 


f B 1 1 1 

B 


[(s + a) s=0j 

(s + a) 

s«0 


(t - t) (2-9) 


which at time t yields 


. «. B T 

y 2 (t) — r • 


(2-10) 


Since at time t the system output is zero, then 

O-VO + yjW-^-V , 


( 2 - 11 ) 


or 


t 

o 


B T 
A 


(2-12) 


which is identical to (2-6). 

A More General Viewpoint . The identical result shown in (2-6) and 
(2-12) is not a coincidence. Examining (2-1) and (2-2), we see that the 
threshold T g drops out from the analysis. That is why T g does not appear 
in the expression for t^. From the physical point of view, corresponding 
to each criterion for bias compensation, there is a threshold level as 
shown in Fig. 2-4. The period of the bias compensation pulses is simply 
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the time required for the ramp due to the bias drift to achieve the magni- 
tude of the system response due to a single torquing pulse. At the instant 
of pulse application, we therefore have 



( 2 - 13 ) 


t 

o 


A 


( 2 - 14 ) 


EFFECT OF WAVE SHAPE UPON SENSOR ACCURACY 

In the previous chapter the bias error, which is generated by the 
electronics network under quiescent conditions, was examined. Under con- 
ditions of variations in attitude the additional error due to the nonideal 
shape of the torquing pulses will now be determined. Since the pulses must 
be calibrated so that the effect of each one upon the sensor will be pre- 
cisely known, the sensitivity of the sensor to variations in the pulse shape 
is a prime consideration. Utilizing the effects of both ideal and realizable 
pulses upon the sensor, a comparison of the result as a function of pulse area 
forms the basis of this sensitivity analysis. 

Property of Dynamic System with Staircase Output Response 

In general, a system containing a pulse width modulator is nonlinear. 
However, it will be shown that, because of the way the output signal is fed- 
back, the dynamics of the strapdown sensor rebalance loop is linear as long 
as the pulse width modulator is not driven beyond its saturation limit. Thus, 
linear analysis techniques can be applied for the study of the rebalance loop. 

Consider a time invariant system G(s) with the associated impulse re- 
sponse g(t). The output y(t) and input u(t) are related by the real convolu- 
tion integral 

t 

•/ 


y(t) 


g(t - t) u(t) dt . 


( 2 - 15 ) 
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Next, attaching a sampler and a zero order hold at the output of 
G(s) results in a new system G'(s) as shown in Fig. 2-7. This new system 
has a staircase shape output response. The impulse response of this system 
is given by 

g’(t) - g(iT) , iT < t < (i + 1)T . (2-16) 


Note that g' (t) -*■ g(t) for t -► o. The general output y(t) responding to an 
input u(t) is 


y'(t) 


f 


g' (t - t) u ( t ) dr . 


(2-17) 


At the discrete time t * nT , n-0 ; 1, 2, -- this now becomes 


y’(t) 


r 


g’(t - t) u ( t ) dx 



g* (nT - t) u (t ) dx . 


But using (2-16), yields 


n 


y’(t) - £ 
^■-00 


IT 

I g’ (nT - i - 1 T) u (x ) dx . 
(i-l)T 


Since g*(nT - i - IT) does not involve the variable x, it can be taken out- 
side of the integral sign, i.e. , 


n 


iT 


y’(t) 


E 


g' (nT - i - IT) 


/ 

“ / (i-.l)T 


u(x) dx 


(2-18) 
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Now define 


A(i) 


f 

(i-l)T 


u(t) dx 


(2-19) 


which is the "input area" over the sampling period ending at t = IT. Then 
(2-18) becomes 

n 

y'(t) = Yj g'< nT " A(i) (2-20) 

j[r=— 


Eq. (2-20) shows that "if the linear system has a staircase output response, 
then the response depends linearly on the input area over each sampling per- 
iod, but not on the input wave shape during that period". 

The Equivalent Rebalance Loop Configuration 

Consider a typical pickoff-electronics block for width-modulated bi- 
nary pulses. The input-output relationship is shown in Fig. 2-8a. This 
relationship can be represented by an equivalent block diagram shown in Fig. 
2-8b. Therefore the overall dynamics of the rebalance loop can be represented 
by a configuration shown in Fig. 2-9 which is equivalent to Fig. 2-1. From 
the conclusion obtained in the last section and expressed in Eq. (2-20), the 
dynamics of this system is linear. 

Effect of Time Constant on Pulse Area 

Since pulse area is the determining factor in the system response, its 
sensitivity for the circuit parameter variations is important. Here we shall 
consider the variation in pulse area caused by a variation of the circuit time 
constant. 

Let an ideal pulse be represented by 

f(t) = u(t) - u(t - -j) (2-21) 

where u(t) is *-he unit step function and the pulse period is T. The area of 


Float 

Angle 

a 


Plckoff & 
Electronics 


Current Pulses 



(a) 





(b) 


Current 

Pulses 


Figure 2-8. The Equivalent Pickoff-Electronics Block Diagrams 





In] 

Ra 



Figure 2-9. The Equivalent 
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this pulse is 


a T 

A O — 
2 


Next, consider a pulse with finite and equal, but non-zero, rise and fall 
times as shown in Fig. 2-10. This pulse is represented by 


£’(t) - 1 - e T ju(t) - 1 - 




where x is the time constant, and its area over one period is given by 


f-[ 


e - e x 


Now define an "area ratio", A.R. , as a measure of area deviation, 

. _ Actual area 

A.R. *= vi — ; 

Ideal area 


Then (2-22) and (2-24) give 


A.R. = 1 - 


(2-26) 


Table 2-2 shows the area ratios for different values of circuit time constants. 

Table 2-2. Area Ratios 


Time Constant x sec. 

Area Ratio 

io“ 3 

0.52269 

1ST* 

0.99866 

5 x 10 -5 

0.99999 

2 x 10" 6 

.0000 

l<f 6 

1.0000 
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In order to form an effective measurement of area loss, the loss 
per unit time is used since this loss will be observed on each input pulse. 


. , , , 2 (1-Area Ratio) 

Area Loss/second = — 5 - *- 


(2-27) 


For example, an Area Ratio of .99999 for a torquing pulse with a width of 
two milliseconds will yield an Area Loss of 0.01 pulse per second. 


, 2(1. - .99999) 

Area Loss/second = — » — - ■ 10 

2 x 10 


(2-28) 


This loss will accumulate to the equivalent loss of one pulse in one hundred 
seconds, assuming full torquing of one polarity over that interval. This 
loss of one pulse for typical strapdown sensor loop resolution will result in 
errors on the order of 1 arc second. 

Effect of Pulse Shape Upon the Output 

Since the effects of rise and fall times upon the input pulse area 
are now known, the effect upon the output of the system is needed to com- 
plete the analysis. Using the system, 

G(s) « - r (2-29) 

s(s + 1(T) 


and the ideal and nonideal pulse inputs which in Laplace notation are rep- 
resented by 


R(s) 



-Ts/2 


1 


9 


(2-30) 


and 


R(s) 


1 

s 


T 

1+TS 



) 


(2-31) 


respectively. The variations in outpu' 4 ue to the nonzero rise and fall 
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times can be determined. 

In the response of the system to the ideal input pulse is given by 


C(s) 


(1 - e~ T ° /2 ) 
[s 2 (s + a)] 


(2-32) 


Taking the inverse Laplace of C(s) yields 


c (t) - (Iff 3 t - l(f 6 + 10" 6 e“ 10 t ) u (t ) 


- [l(f 3 (t - T/2) - 10' 6 


+ 10" 6 e “ 103(t - T/2) ] «(t - T/2) 


(2-33) 


Likewise, the system output for a nonideal pulse input is 


C (s) 


(1 - e" Ts/2 ) _ (1 - e ~ Ts/2 ) 

s 2 (s + 10 3 ) s(s + g)(s + 10 3 ) 


(2-34) 


which in the time domain is 


[■ 


C(t) * 1 10" 3 t - 10“ 6 + 10' 6 e" 10 1 - X 10 -3 


3 i 

in -3 -10 J t t 1 

10 r e T e 

13 + 3 u (t) 

1 - t 10 J tlO -1 1 


l<f 6 - 10- 3 (t - b - io- 6 c* 10J<t - 7> 


+ TlO 


“10 J (t - j) 

-3 10 t e 

1 - TlO 3 


t e 


- i(t - f) 


TiO J - 1 


u (t - 


I) 

2 J 


(2-35) 
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The responses for a system when the nonideal pulse area very closely 
approximates the ideal pulse area are shown in Table 2-3 and Figures 2-11 
and 2-12. 

Also, the responses of the system are shown in Table 2-4 and Figures 
2-33 and 2-14 for an input pulse of an area which is about the same as the 
Ideal area, but the rise and fall times are far from zero. 

Remarks 

From the preceding analysis and examples, the system output can be 
seen to be insensitive to variations in pulse shape provided the pulse area 
is constant. Also, the effects of the individual area deviations seem small 
relative to the overall responpe, but the errors due to the area differences 
occur with each pulse causing a cumulative effect which can be seen to be 
quite large for small individual pulse area errors. For this reason, the 
pulse generator must be capable of creating pulses with highly constant areas 
under changing conditions. Since this error is cumulative, as is the bias 
error, it also must be corrected by updating the computer at regular inter- 
vals. This, however, presents no new problem since other errors also neces- 
sitate the correction. 


CONCLUSIONS 

The results of the analysis presented in this chapter are useful in 
three different ways: (1) they allow an estimation of system errors for 
specific circuit designs and configurations, (2) they provide a criterion 
of allowable dc bias and waveshape variation for the circuit designer, and 
(3) when the sensor error is found by recalibration, they provide a means of 
determining how the error can be compensated by the digital computer. 

Additional investigation is needed in the areas of dynamic error com- 
pensation and reliability of the ■'ndividual functional units in the sensor 
loop. The reliability is of prime importance due to the system complexity. 
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TABLE 2-3 

SYSTEM RESPONSE UNDER IDEAL AND NEAR-IDEAL PULSING CONDITIONS 


Ideal 

Input 

Actual 

Input 

Time 

Ideal 

Output 

Actual 

Output 

1.0 

0.0 

0.0 

.3678 

X 

10" 6 

.3082 

x 10" 6 

1.0 

0.99996 

0.001 

1.1353 

X 

10- 6 

1.108 

x 10" 6 

1.0 

0.99999 

0.002 

2.0497 

X 

10- 6 

1.955 

x 10~° 

1.0 

1.0 

0.003 

3.018 

X 

ID' 6 

2.9203 

x 10" 6 

1.0 

1.0 

0.004 

4.0067 

X 

l<f 6 

3.9074 

x 10 -6 

1.0 

1.0 

0.005 

5.0024 

X 

lO' 6 

4.9026 

x 10~ 6 

1.0 

1.0 

0.006 

6.0000 

X 

10~ 6 

5.9000 

x 10" 6 

1.0 

1.0 

0.007 

7.0000 

X 

10- 6 

6.9000 

x 10~ 6 

1.0 

1.0 

0.008 

8.0000 

X 

lO' 6 

7.9000 

x 10“ 6 

1.0 

1.0 

0.009 

9.0000 

X 

10- 6 

8.9000 

x 10" 6 

0.0 

1.0 

0.010 

9.6321 

X 

10- 6 

9.6264 

x 10~ 6 

0.0 

0.00004 

0.011 

9.86467 

X 

10- 6 

9.8621 

x 10~ 6 

0.0 

0.00001 

0.012 

9.9502 

X 

lO -6 

9.9500 

x 10 -6 

0.0 

0.0 

0.013 

9.98i68 

X 

10" 6 

9.98160 

x 10 -6 


, n -4 

t ■ 10 sec. 

-2 

x * 2 x 10 sec. 





0 0.01 


Time 


Figure 2-11. Ideal and Near-Ideal Area Conservation Input Pulses. 



Time 


Figure 2-12. System Response Under Ideal and Near-Ideal Area Conservation Input Conditions. 
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TABLE 2-4 

SYSTEM RESPONSE UNDER IDEAL AND POOR PULSING CONDITIONS 


Ideal Actual Time Ideal Actual 

Input Input Output Output 


1.0 

0.0 

0.0 

0.0 



0.0 



1.0 

0.63212 

0.0001 

0.048 

X 

l < f 6 

0.00132 

X 

ICf 6 

1-0 

0.86467 

0.0002 

0.187 

X 

,< f 6 

0.0603 

X 

Kf 6 

1.0 

0.95022 

0.0003 

0.408 

X 

10- 6 

0.1489 

X 

icf 6 

1.0 

0.98168 

0.0004 

0.703 

X 

10' 6 

0.3735 

X 

io -6 

1.0 

0.99326 

0.0005 

1.065 

X 

l ( f 6 

0.6737 

X 

io -6 

1.0 

0.99752 

0.0006 

1.488 

X 

10- 6 

1.036 

X 

l ( f 6 

1.0 

0.99908 

0.0007 

1.966 

X 

10- 6 

1.462 

X 

io -6 

1.0 

0.99966 

0.0008 

2.994 

X 

l < f 6 

1.942 

X 

io - 6 

1.0 

0.99992 

0.0009 

3.065 

X 

10- 6 

2.472 

X 

io - 6 

0.0 

0.99999 

0.0010 

3.678 

X 

10- 6 

3.046 

X 

io *- 6 

0.0 

0.36788 

0.0011 

4.280 

X 

ID ' 6 

3.662 

X 

io ‘ 6 

0.0 

0.13533 

0.0012 

4.825 

X 

icf 6 

4.227 

X 

io - 6 

0.0 

0.04978 

0.0013 

5.457 

X 

l ( f 6 

4.799 

X 

io - 6 

0.0 

0.018316 

0.0014 

5.763 

X 

icf 6 

5.287 

X 

io - 6 

0.0 

0,006738 

0.0015 

6.166 

X 

10~ 6 

5.712 

X 

io -6 

0.0 

0.000576 

0.0016 

6.531 

X 

icf 6 

6.130 

X 

icf 6 

0.0 

0.0 

0.0017 

6.862 

X 

10- 6 

6.520 

X 

io " 6 

0.0 

0.0 

0.0018 

7.158 

X 

10- 6 

6.861 

X 

io " 6 

0.0 

0.0 

0.0019 

7.430 

X 

10- 6 

7.172 

X 

icf 6 

0.0 

0.0 

0.P02C 

7.675 

X 

10- 6 

7.442 

X 

io " 6 

0.0 

0.0 

O.CC?, 

7.895 

X 

l<f 6 

7.805 

X 

io - 6 



Figure 2-13. Ideal and Poor Area Conservation Input Pulses. 




Figure 2-14. System Response Under Ideal and Poor Area Conservation Input Conditions. 
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CHAPTER 3 


STRAPDOWN PLATFORMS USING REDUNDANT 


TWO-DEGREE-OF-FREEDOM GYROS 


J.C. Hung and LeSi Han 


' N72-?7676 


ABSTRACT 


This chapter presents a new concept of high 
reliability strapdown attitude sensing systems for 
space vehicles. Each system utilizes a set of re- 
dundant two-degree-of-freedom gyros. An optimum 
system configuration is obtained for maximum system 
reliability and the best measurement accuracy. Im- 
proved accuracy of the final data is obtained by 
using the least-square data reduction technique. 

Each system possesses a "sensor performance manage- 
ment" feature which is capable of failure detection, 
faulty gyro identification, system reconfiguration, 
and possibly, sensor recalibration. Improvement in 
reliability as compared to other types of strapdown 
systems is demonstrated. Details of the development 
are described in terms of a system containing four 
gyros. Results from systems containing three, five, 
and six gyros are included in the Appendix of this 
chapter. 


INTRODUCTION 

This chapter proposes a new redundancy scheme for high reliability 
strapdown navigation systems. The scheme involves the use of two-degree-of- 
rreedom (TDF) gyros rather than SDF gyros for attitude sensing. The advan- 
tages of using TDF gyros as compared to SDF gyros are 1) higher reliability 
for the same number of gyros used, 2) less effect of sensor inaccuracy, 3) 
simpler computation for data reduction, and 4) provision of more redundant 
measurements for better data reduction. In this chapter, the system reli- 
ability of the proposed scheme will be presented; comparison will be made 
with respect to other known schemes; a discussion on the optimum configura- 
tions for TDF redundant sensors will be given; and the methods of failure 
detection, faulty sensor identification, and system reconfiguration will be 
described in detail. 
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SYSTEM RELIABILITY 

Each TDF gyro senses attitude motion in two orthogonal directions. 
Therefore an attitude sensing system composed of n TDF gyros will operate 
satisfactorily if only 2 out of n gyros are in good condition as compared to 
3 out of n for the system using SDF gyros. Thus, with the same number of gyros, 
the system using TDF gyros provides more redundancy and thus higher reliability 
than the one using SDF gyros. 

Analytically, reliability is defined as the probability of success. 

Assume all individual gyros have equal reliability 

r «= e" Xt (3-1) 

where A is a positive real number and t represents time. Eq. (3-1) shows that 
reliability is a monotonically decreasing function of time. For an attitude 
sensing system formed by n TDF gyros the system reliability is given by 

R = ( n ) r 11 + ( n 1 )r n " 1 (l - r) + + (") r 2 (1 - r) n " 2 (3-2) 

n n n-1 i 

a similar system formed by n SDF gyros has a system reliability of 

R’ = ( n ) r 11 + ( n .) r 11 " 1 (1 - r) + + (*) r 3 (1 - r) n " 3 

n n n-1 -3 

“ R - (") r 2 (1 - r) n ' 2 (3-3) 

n z 

A plot of R and R* , for n = 2 to 6, versus the normalized time T = At is 
r n n 

shown in Fig. 3-1. It is evident that higher reliability is associated with 
TDF gyro strapdown systems. Note that in the above discussion it has been 
assumed that the reliability r of each individual SDF gyro is the same as 
that of each TDF gyro. This assumption is reasonable from a hardware con- 
struction point of view. 

OPTIMUM SYSTEM CONFIGURATIONS 

Three considerations are involved in arranging gyros for optimum system 
coiif igurations. The first and the most important consideration is the maximum 
system reliability. For example, a group of four TDF gyros A, B, C, and D can 
be oriented with their sensitive axes directed in a way symbolically shown in 
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Fig. 3-2. In this figure, the circles .labeled A1 and A2 on x and y axes, re- 
spectively, represent the directions of two sensitive axes of gyro A. Similar 
representation is used for all other gyros. This arrangement is not optimum, 
since if gyros A and C fail the remaining gyros B and D cannot provide a full 
attitude sensing. 

The second consideration is the effect of gyro orientation on measure- 
ment error introduced by sensor inaccuracy, kcferring to Fig. 3-3, let u be 
the desired vector quantity and S is the direction of sensor's sensitive axis 
which makes an angle a from vector w. The measured quantity is 

m ■ w cosoi + e 

where e represents the measurement error. The "normalized measurement error" 
is given by 

e = — (3-4) 

N u cosa 

Eq. (3-4) shows that this error is the least when a = 0. In other words, if we 
want to measure the vehicle attitude motion along three vehicle coordinates, 
the measurement error is the least if all sensitive axes of sensors are directed 
along vehicle coordinates. 

The last consideration is the simplicity of computation. Simpler com- 
putation requirements result in less computation errors and shorter computation 
time. In a trapdown navigation system usually the attitude motion along 
vehicle axes is measured . Therefore a TDF gyro whose sensitive axes are along 
the vehicle axes senses directly the attitude motion along those two vehicle 
axes, no computation is needed. 

It is assumed that if any one or both outputs of a gyro fail, the 
entire gyro is taken out of service. 

Guided by the three considerations and the assumption above, optimum 
system configuration for high reliability strapdown systems using three, 
four, five and six TDF gyros are given below, where the desired quantities 
are vehicle rates w x , oj^, and u> z along three vehicle axes. The arrangements 
are also symbolically shown in Fig. 3-4. 
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Gyro Directions of sensitive axes 
A x and y axes 

B y and z axes 

2 z and x axes 

D y and z axes 


Figure 3-2. A Non-Optimum System Configuration. 
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S 



Desired signal: cu 

Measured signal: od cos a + e 

e 

Normalized measurement error: ~' " co8 g 


Figure 3-3. Normalized Measurement Error. 




Figure 3-4. Optimum System Configuration. 
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3 TDF Gyro8 

The arrangement and measurements for a 3 gyro system are shown in 
Table 3-1. The measurement equation is 

m(3) - R < 3) w , (3-5) 


where 


co 


CO 


(0 


CO 


(3-6) 


and 


m(3) - [nj m 2 m^ m 5 m 6 ] 


R(3) 


1 0 0 0 0 1 1 
0 110 0 0 
0 0 0 1 1 0 


(3-7) 


The symbol [ ] denotes the transpose of a matrix. 

4 TDF Gyros 

The arrangement and measurements for a system composed of four gyros 
are shown in Table 3-2. The measurement equation is 


m(4) ■ R(4) co 


(3-8) 


where 


m(4) - [m^ m 2 m g m^ m^ m g m ? m g ] T 

1000011 0** T 
R(4) - 0110000c 

0001100 cj 

Throughout this paper c « cos45° - 0.707. 


(3-9) 


5 


Direction of 
sensitive axis 


Measurement 


Al, 

x 

direction 

CM 

< 

y 

direction 

Bl, 

y 

direction 

B2, 

z 

direction 

o 

h-* 

z 

direction 

C2, 

X 

direction 


m, *= o 
1 x 


m. *» cd 
4 z 


m, *= cd 
6 x 


Table 3-1. Arrangement for 3-TDF-gyro system 
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5 TDF Gyros 

The arrangement and measurements for a 5 gyro system are shown in 
Table 3-3. The measurement equation is 


m(5) c R(5) a) 


(3-10) 


where 


m(5) = [n^ m 2 m 3 m^ m 5 m fi m y m g m g m 1() ] 


R(5) 


100001100c 

OllOOOOclO 

OOOllOOcOc 


1 T 


(3-11) 


6 TDF Gyros 

The arrangement and measurements for a 6 gyro system are shown in 
Table 3-4. The measurement equation is 


m(6) - R(6) w 


(3-12) 


where 


m(6) = [m^ m 2 m g 


m 12 ] 


lOOOOllOOcOc 
R(5) = OllOOOOclOOc 
OOOllOOcOclOj 


i T 


(3-13) 


Notice that in any of these configurations, any two good gyros can provide a 
complete attitude sensing. 


SENSOR PERFORMANCE MANAGEMENT 

A high reliability navigation system employing a redundancy concept 
must be capable of "failure detection", "faulty sensor identification , 
"system reconfiguration", and possibly "sensor recalibration". All these 
functions can be handled by a centralized control called "sensor performance 


Gyro 


Direction of 
sensitive axis 


Measurement 





A 

Al, x direction 

m, = cd 

1 x 

A2, y direction 

a 

N> 

ii 

B 

Bl, y direction 

m 0 = cd 

3 y 

B2, z direction 

a 

•p* 

ii 

e 

N 

C 

Cl, z direction 

m_ = cd 
5 z 

C2, x direction 

m, = cd 
6 x 

D 

Dl, x direction 

B 

II 

D2, y = z direction 

m 0 = 0.707 cd + 0.707cd 
8 y z 

E 

El, y direction 

m_ = cd 
9 y 

E2, z = x direction 

m, = 0*707 cd + 0.707cd 
10 z x 
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management" (SPM) . This management task is performed by a computer software. 
Two SPM schemes are shown in Figures 3-5 and 3-6 in the form of operational 
flow diagrams. The first scheme involves failure detection, faulty sensor 
identification, and system reconfiguration. The second scheme involves the 
additional ability of determining whether the failure is a soft failure and 
the ability of sensor recalibration in case of a soft failure. 

When the system consists of n gyros, the total number of states is 
2 n , and 2 n - (”) - (”) of them allow faulty gyro identification. 

In the following sections, details of the SPM concept are developed for 
a system containing four TDF gyros. Details for other systems are included 
in the Appendix of this chapter. The development is described in three sepa- 
rate parts, namely, failure detection and faulty gyro identification, system 
reconfiguration and data reduction, and sensor recalibration. 

DETECTION AND IDENTIFICATION 

The function of failure detection is to find out if a failure has oc- 
curred somewhere in the system while the function of faulty gyro identification 
is to identify the faulty gyro. When the system contains four TDF gyros, there 
is a total of sixteen possible states including the cases when all gyros are 
good, one is bad, two are bad, etc. It will be shown that although all states 
allow failure detection, not all of them allow faulty gyro identification. 

Parity Check Equations 

When all four TDF gyros are normal, the following conditions exist as 
can be seen from the measurement equation (3-8). 

i) m 6 - m 7 “ 0 

ii) m g - c m 3 - c m^ = 0 

iii) m^ - m^ “ 0 



m 2 - m 3 = 0 


iv) 

v) 

vi) 


✓ 


(3-14) 
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The "approximately equal" sign accounts for the normal measurement error. 
Assume that no two faulty gyros have identical faults; or, if they do have 
identical faults, assume that both do not become bad at the same time. Then 
one or more of the six equations above will not be satisfied if any failure 
occurs. These six equations are called "parity check equations". 

Truth Table for Faulty Gyro Identification 

Let us associate with the ith parity check equation a binary valued 
quantity 1^. ¥L ± equals "0" when the equation is satisfied and equals "1" 
when not satisfied. 

Assume that when a TDF gyro fails, measurements from both axes are 
erroneous. The validity of this assumption is subject to discussion and 
further exploration. Now we are in a position to construct a truth table 
for faulty gyro identification. Suppose that gyros A and are bad, then 
measurements m^, m 2 , m^ and m^ are not consistent. Under this condition 
all parity check values are equal to 1, except k^ which is equal to 0. 
Reasoning in this manner we obtain a truth table which is shown in Table 
3-5. 

Examining the table shows that there are twelve distinguishable system 
states. Eleven of them, which include the events of not more than two failed 
gyros, allow both failure detection and faulty gyro identification. The last 
state, which includes the group of five events where more than two gyros are 
failed, allows failure detection only but not the identification of faulty 
gyros. 


DATA REDUCTION AND SYSTEM RECONFIGURATION 

Aside from increasing the system reliability v:ne use of redundant gyros 
provides redundant measurements which can be used by means of a data reduction 
technique to reduce the effect of errors associated with individual gyros. 

The technique used will be the "least-square data reduction". Since the 
system reconfiguration is accomplished by means of software control, it can 
readily be incorporated into the data reduction program. 

Data Reduction 

Rewrite Eq. (3-8) in its explicit form and omit the parentheses asso- 
ciated with ffi and R for simplicity. 
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Table 3-5. Truth fable for failure detection and faulty 
gyro identification, with reconfiguration 
solution program indicated, for 4-gyro system 
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(3-15) 


In this equation the measurement m is available from sensor outputs while 
the attitude rate u is desired. Since R has more rows than columns, no 
simple to can be found to satisfy this equation. Instead, any selected _to 
will result in a solution error 

£ “ R to - m (3-16) • 

The least-square error criterion is to choose a jo such that the quality 

I ■ £*£ * (Rxo - m)T (Rto - m) (3-17) 


is minimized. Taking the gradient of (3-17) with respect to jo and setting 
the result to zero gives 

VI - R T (Pxo - m) - 0 (3-18) 


Solving for u yields 
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0 0 it 
u 12 12 12 12 y 
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0 12 12 12 12 0 0 f 


(PO) 


which is the desired data reduction equation readily be written into a 
computer program. Let us label this equation and its associated program PO. 

It should be mentioned that if the statistics of the sensor inaccuracies 
are known, more sophisticated techniques can be used for better results. 


System Reconfiguratio n 


The system reconfiguration need not be done by hardware modification. 
Instead, it is conveniently and rapidly done by changing the software program, 
which can be executed by the computer. Consider the case that gyro A is found 
faulty. Under this condition measurements m^ and m 2 will not be used in the 
data reduction. This is done simply by eliminating m^ and m 2 rows from Eq. 
(3-14), giving 
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(3-19) 


Denoting the reduced measurement vector m^ and measurement matrix R^, the 
new least-square data reduction equation should be 
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- = < r i T r i ) ” 1 r i T -1 * (P1) 

To generalize, the system reconfigurations for the first ten failure 
conditions can be accomplished by using the appropriate measurement vectors 
and measurement matrices in the data reduction equation 

- " (R i T V** 1 R i T m i » < pi > 

i ■ 1 to 10. For the last five failure conditions, no system reconfiguration 
is possible, and the computer should report a total-system failure. Table 3-6 
lists all the measurement vectors m^ and measurement matrices for various 
failure conditions. 

It should be mentioned that there are other possible data reduction 
techniques. For example, if the statistics of sensor inaccuracies are known, 
more sophisticated techniques can be employed for better data reduction. 

SENSOR RECALIBRATION 

If a faulty gyro has a soft failure, such as a bias error, a change 
of scale factor, cr showing erratic readings about a mean value, it can be 
recalibrated. The faulty gyro is taken out of service by system reconfigu- 
ration and put into a fault determination mode. During this mode, outputs 
of this faulty gyro are compared with the outputs of good gyros to determine 
the type of fault which has occurred and the amount of correction needed. The 
correction is made by changing parameters in the software. The system is then 
reconfigurated once. 

While the functions of failure detection, faulty gyro identification, 
and systc-u. reconfiguration are indispensable in sensor performance management, 
the function of sensor recalibration is not. The inclusion of the latter 
maximizes the system reliability at the expense of an increased computer re- 
quirement, which may be a worthwhile investment. 
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Table 3-6. Vectors and matrices for reconfiguration 
solution program for 4-gyro system. (c=0.707) 
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APPENDIX 


SPM DETAILS FOR SYSTEMS HAVING 
THREE, FIVE, AND SIX TDF GYROS 


3-Gyro System 

The set of three parity check equations for the 3-gyro system is 
given by 



(3-20) 


The truth table for the identification of faulty gyros are shown in Table 
3-7 with reconfiguration solution programs indicated. The data reduction 
equation for the normal condition, where all gyros are good, is given by 
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Table 3-7. Truth table for detection and faulty gyro identification 
with reconfiguration solution program indicated, for the 
3-gyro system. 
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where m = [m^, m^ m^, m,., m^,] and 


10 0 
0 10 
0 10 
0 0 1 
0 0 1 
10 0 

There are seven abnormal conditions, of which only three are 
identifiable. Each reconfiguration solution program for data reduction 
has the form 


T -IT 

w = (R t \) R ± Hi 


1, 2, 3 , 


(PI) 


where the measurement vectors m^ and measurement matrices are listed 
in Table 3-8. 

5-Gyro System 

The set of ten parity check equations for the 5-gyro system is 
given by 

\ 

m g + c m 7 - m 1Q - c m ? = 0 

■lo ' c V c *6 “ 0 

m 6 ' m 7 * 0 
m g - m^ - 0 
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8 3 4 
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4 5 
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(3-21) 
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Table 3-*8. Vectors and matrices for reconfiguration solution 
programs for 3-gyro system. 
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The truth table for the identification of faulty gyros are shown in Table 
3-9, where the reconfiguration solution program for each abnormal condition 
is also shown. The data reduction equation for the normal condition where 
all gyros are good again has the form of Eq. (PO). But here 
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(3-22) 


There are thirty-one abnormal conditions, of which only twenty-five 
are identifiable. Each reconfiguration solution program for data reduction 
has the form 


6) 




l i =1 


i - 1 to 25 


(Pi) 


where the measurement vectors m. and measurement matrices R. are listed in 

—1 l 

Table 3-10. 


6-Gyro System 

The set of fifteen parity check equations for the 6-gyro system is 
given by 
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Table 3-9. Truth table for failure detection and faulty gyro 

identification, with reconfiguration solution programs 
indicated for 5-gyro system. 
























































Reconf igurat ion 
Solution Program 


m.f and are obtained from m and R by 
removing the following rows from Eq. (3-22) 
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Table 3-10. Vector m. and matrices R^ for reconfiguration programs 
for the 5-gyro system. 
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The truth table for the identification of faulty gyros are shown in Table 
3-11 where the reconfiguration solution program for each abnormal condition 
is also shown. The data reduction equation for the normal condition where all 
gyros are good is again given by Eq. (P0). But here 
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(3-24) 


There are sixty-three abnormal conditions, of which only fifty-six 
are identifiable. Each reconfiguration solution program for data reduction 
has the form 

u - (R. T R,)" 1 R. T m, i - 1- 56 (Pi) 

— i i i — i 

where the measurement vectors m, and measurement matrices R, are listed in 
Table 3-12. 


Recon f Igoration 
Solution 
Frogran 


tnj Rj are obtained from * and It by 
removing the following rows from Eq. 0-21) 


F4 

?. e 

F5 

9, 10 

F6 

11, 12 


1 . 2 , 7 . 8 


1, 2. 9. 10 


1 . 2 . 11 , 12 


3, *. 5, 6 

3. 4. 

2. 

8 

3. 4, 

9, 

10 

3, 4, 

n. 

12 

5, 6, 

7. 

8 

5, 6. 

9, 

10 

5, 6, 

11, 

12 

7. 8, 


10 

7, 8, 

11. 

12 

9, 10 

. U 

, 12 


1, 2. 3. 4, 7. 8 


1, 2. 3, 4. 11. 12 


1, 2, 5. 6. 7. 8 


1. 2, 5, 6. 

9. 

10 

1, 2. 5. 6, 

11. 

12 

1. 2. 7. 8. 

s, 

10 

1. 2, 7. 8, 

11, 

12 

. ... . 1, 2, 9, 10 

i 11 

. 12 

3. 4. 5. 6. 

7, 

8 

3, 4, 5, 6, 

9, 

10 

3. 4, 5. 6, 

11. 

12 

3, 4, 7, 8, 

9, 

10 

3, 4, 7, 8, 

11, 

12 

3, 4, 9, 10, 

, 11 

. 12 

5, 6. 7, 8, 

9, 

10 

J, 6. 7. 8. 

11. 

12 


7, 8, 9, 10, 11. 12 


1. 2. 3. 4. 3. 6. 9. 10 


1. 2. 3. 4. 7. 8. 1 


1. 2, 3, 4. 9. 10. 1 


1, 2, 5. 6. 7, 8. 9, 10 


1, 2, 5, 6, 7, 8, 11. 12 


1, 2, S, 6, 9, 10, 11, 12 


1, 2, 7, 8. 9, 10, 11, 12 


3. 4, 3, 

_ 6 . 

_2. 

3. 9, 10 

3. 4, S. 

6. 

7, 

8, 11. 12 

3. 4, 5, 

6. 9. 

10, 11, 12 

3. 4. 7. 

8, 

9. 

10, 11, 12 

5. 6, 7, 

3, 

9. 

10. 11, 12 


Tula 3-12. Vectora «j and aatricaa *j for reconfiguration 
program* for the 6-gyro aystem. 
















































CHAPTER 4 


DDA REALIZATIONS OF ATTITUDE ALGORITHMS 
J.C. Hung and Y. Yeh 
ABSTRACT 



This chapter presents 1) a review of six atti- 
tude transformation algorithms, and 2) the realiza- 
tions of these algorithms by digital differential 
analyzers (DDA’s). The number of DDA's for each 
algorithm is either smaller than, or equal to that 
obtained by others. 


INTRODUCTION 

The use of strapdown platforms for the guidance and control of 
space vehicles has received a good deal of attention in recent years. In 
using a strapdown platform, the attitude of a space vehicle is measured 
along a coordinate frame fixed to the vehicle body. However, it is commonly 
desired to determine the attitude of the space vehicle with respect to a 
different coordinate frame, such as an inertial coordinate frame or an earth 
fixed coordinate frame. Therefore, a transformation of attitude information 
from the vehicle body frame to a desired coordinate frame is needed. Several 
different attitude transformation algorithms are available. The most common 
ones are (1) direction cosine algorithm, (2) Euler angle algorithm, (3) Euler's 
four parameter algorithms, (4) quaternion algorithm, (5) Cayley-Klein algori- 
thm, and (6) Vector representation algorithm. 

An analog computer, a digital computer, or a group of digital differ- 
ential analyzers (DDA's) can be used to implement the attitude algorithms. 

The advantages ot using DDA's include high accuracy and compact size. 

This chapter begins with a review of the analytics for the six atti- 
tude algorithms. Then, each algorithm is realized by a group of digital 
differential analyzers (DDA’s). It will be shown that the number of DDA's 

required for the realization of each algorithm is either smaller than, or 

1 2 

equal to that obtained by others. * The different numbers of required DDA's 
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for various algorithms, also serve as an indirect comparison of computation 
requirements if the algorithms are to be lealized by digital computer. 

THE DDA 

DDA is a device which has two incremental inputs Ax^ and Ay^, one 
memory memorizing the quantity and one output Az^. The device per- 

forms the following mathematical operation 

Az k " <y k-i + iy k ) lx k (4_1) 

The following symbol is used to represent a DDA. By connecting a number of 



Figure 4-1. The DDA 


DDA's together in conjunction with summing devices, various computation algo- 
rithms can be realized. 


DIRECTION COSINE ALGORITHMS 
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Consider two orthogonal coordinate systems S^ and S 2 in a three dimen- 
sional space: 
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For any vector V in the space, it can be expressed in terms of coordinate 

3 

components in "vector-analysis" form 


V B i v^ + j Vj + k v^ e i'v| + j'v^ + k'v^ 


(4-2) 


By the property of dot product v^ can be expressed in terms of Vj as 
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In "vector-matrix" form^, we have 
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where v is a column vector whose three elements are components of V along co- 
ordinate system and v’ is a column vector whose three elements are components 
of V along coordinate system S 2 The matrix C is called the "direction cosine 
matrix", whose elements are cosines of angles between two sets of unit vectors 
of and S 2 . Thus (4-4) can be written as 


C v 


(4-5) 


A similar derivation by expressing v^ in terms of v^ gives 


C T v' 


(4-6) 



Using both (4-5) and (4-6), 


T 

v « C C v 


that is 


T 

C C - I 


(4- 


Hence, C is an orthogonal matrix. 
(4— T' ) says that 


3 

I 

k-1 


C mk C kn 


Denote the elements in C by C 


i y 


Eq 


1 if m " n 
0 if m ^ n 


(4* 


which is a set of six constraint equations. Therefore, among the nine 
elements of C, only three are independent. 

An Example: A two dimensional case. 
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t < 


» A 

A 

A 

A 



v 'l 

m 

i 1 • 

i 

i* * 

j' 


V 1 

• V 2 


A 

[ y • 

A 

i 

A 

y • 

A 

j _ 


V 2 



COS0 


cos (90+0)1 

V 1 


*= 









_ cos (90+0) 

COS0 


■ 

V 2 



COS0 


sin0 

• 


V 1 


B 









-sin0 


cos0 

. 


^2 


Relationship Between Two Sets of Unit Vectors 

Consider a vector V making angles a, B, and y with the set of unit 

A A A 

vectors i, j, and k, respectively, as represented by 

AAA 

V = i cosa + j cosB + k cosy (4-9) 



Figure 4-3. Angles a, B, and y. 

2 2 2 

Since cos a + cos B + cos y = 1, V is also a unit vector. If we let V = i f , 

A' A 

j*, and then k' , we have the relationship between the two sets of unit vectors. 
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A A A A A A A A A A » 

i' - (i' • i) i + (i* * j) j + (i' * k) k 

A A A A A A A A A A 

j’ = (j’ • i) i + (j' * j) j + (j’ * k) k > 

A A A A A A A A A A 

k' ■ (k 1 • i) i + (k * * j) j + (k ' • k) k , 


(4-10) 


Representing the two sets of unit vectors by 



< 


A 

i’ 

u - 

A 

j 

and U ' “ 

A 

j’ 


■ k ■ 


A 

. k\ 


(4-11) 


(4-10) can be expressed as 

U’ - CU (4-12) 

where C is defined in (4-4). Notice that U and U' are not vectors; their 
elements are unit vectors! 

Now consider any vector V. From (4-2) 



i 


A « 

i’ 

V * (v l v 2 v 3> 

3 

= [v’ v^ vj] 

W A 

r 

T 

V 

. k . 

- 

V’ T 

. k’. 


U 


U' 

or simply 

T 

V « v U « 

v ,T U* 


Using (4-12) 

v ,T U’ - v T 

T T 

U * v C U' 


indicating 

,T 

v ■ 

T_T 
v C 


therefore 

v' - 

C v 



(4-13) 


(4-14) 


which is the same as (4-5). 
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Differential Equation of Direction Cosine Matrix 

A A A 

Let S frame with unit vectors i, j, and k be stationary, and S„ frame 

J* A A A 4 

with unit vectors i', j', and k' be rotating. Let 


u 


CO 


CO 


co 


be the angular velocity of with respect to expressed along S£ axes, and 
C be the "direct cosine matrix" transforming coordinates to coordinates, 
l.e., 


C, 

t 

c . 


A A 


11 

C 12 

C 13 


i • i' 

A A 

i * j' 

A A 

i • k’ 

A A 


21 

C 22 

C 23 


j • i’ 

j * j’ 

A A 

j • k' 

A A 

(4-15) 

31 

C 32 

C 33 * 


, k • i' 

k • y 

k • k'. 



Note that this C is different from that in (4-4) ! One is the transpose of 
the other. 


'11 


'21 


: di' 

= 1 • dT 


« i . dil c 

dt C 13 

; dk' 
dt 

*= i • (j’co - k'co ) 
J z y 


= i • (k’ co - i'co ) 
x z 

= i • (i'co, 

" C 12 U z " C 13“y 


“ C 13“x - C ll“z 

“ C. .co - 
11 y 

. : . di' 

J dt 

c 22 

. i . ill • 

2 dt °23 

. j . dk* 
J dt 


j • (j'co z - k’oo y ) 
C 22“z - C 23“y 


A A 


J * (k'u - I'm ) 

J X z 


C 23“x ’ C 2l“* 


i - 3'» x ) 

C 2l“y ' C 22“x 


• " di 1 

C ■ Ir • — 

31 dt 


'32 


C . ii' 


dt 


c - i . 

33 dt 


k • (j ’w - k'co ) 
J z y 

C- 0 gj - 
32 z 33 y 


k • (k ' co - i ' co ) 
x z 


C 0 ,co - C 01 oo 
33 x 31 z 


k -(i'cOy - j*u x ) 

C 3l“y ' C 32 W x 
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Let 


11 

n 

“12 

C 13 

21 

C 22 

C 23 

31 

C 32 

C 33 ' 


’ 0 

—to 

z 

Q - 

U) 

z 

0 


-0) 

y 

to 

X 


—to 


-w U) 
z y 


0 -w 


to 0 

x 


(4-16) 


be the matrix characterizing the angular velocity of the S£ frame with respect 
to the frame. Then we have 

C - C 8 (4-17) 


Any vector in coordinates can then be transformed to coordinates via 


2.1 


C ^s2 


(4-18) 


In a strapdown navigation system, to , to , and oj are measured by a 

x y z 

?et of strapdown rate gyros. The transformation from v^ to is represented 
by the diagram shown in Fig. 4-4. 



n 

from rate gyros 


Figure 4-4. Direction Cosine Transformation. 
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The symbol 



C 


B 


indicates the matrix multiplication C *■ AB. 


EULER ANGLE TRANSFORMATION ALGORITHM 


Euler Angle Transformation Matrix* * 

Two coordinate systems (X, Y, Z) and (x, y, z) can be related by 
three rotations of coordinates as shown in Fig. 4-5. 


First, rotate along Z axis by an angle $• 

Next, rotate along x^ axis by an angle 6 . 

Finally, rotate along Z2 axis by an angle ip. 

The angle 4), 0, and are called "Euler angles", and the three successive 
transformations are represented by the following equations 


r > 


‘ X 1 


cos<p sin<J> 0 


X 

y l 

- 

-sin<f> cos4> 0 


Y 

. Z 1 


0 0 1 . 
v v ✓ 


z . 

V 


0 

0 

H 


' V 

y 2 

m 

0 cos0 sin9 


y l 

4 


0 -sin8 cos0 

' V 

f 

■ Z 1 


(4-19) 


(4-20) 


79 


X ' 


cost|> 

sin’l' 

0 ‘ 


X 2 

y 

=7 

-sin^ 

COSlJl 

0 


y 2 

z 


0 

0 

1 


Z 2 . 



> 





Thus, the overall transformation is given by 


> « 
X 


' X 

y 

- r 

Y 

z 

« 


. Z . 


( 4 - 21 ) 


( 4 - 22 ) 


where the transformation matrix is 


r - r. r. r. 

ij« 6 <)> 


cosi^ cos^i - cos0 sirnfi sin^ 
-sini{» cos4» - cos0 sin<f> cosi^ 
L sin0 sin<J> 


cosip sin$ + cos0 cos$ sin^ sini^ sin0 

— sin^ sin<J> + cos© cos<p cos^ cosi{> sin0 

-sin0 cosij) cos0 


( 4 - 23 ) 


Euler Angle Differential Equations 

If the xyz-frame is rotating with respect to the XYZ-frame, then the 

Euler angles $, e, and y change with time. The rotation can be described as 

• • • 

a function of their time rates, namely, <f>, 0, and ip. 
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Let u be a vector representing the angular velocity o f the xyz- 
frame with respect to the XYZ-frame. Then 




4 

0 


• 

e 


' o * 

« r r 

Ve 

r— 

-©*• O 

+ r , 
* 

0 

. o - 

+ 

0 

. i . 



sinij; sin0 cost); 0 


♦ 

«r 

cos<J> sin0 -sin^ 0 


• 

0 


cos0 0 1 . 


. * . 


(4-24) 


The time rates of Euler angles can be obtained by inverting (4-24) 

0 
0 

sin0 


* 


sinip 

cos4» 

• 

6 

_ i 

sin0 

zosip sin0 

-sinifi sin0 

.i. 


-sin4» cos0 

• 

-COS0 COSt|) 

• 

t 



T(e) 



’ CO 


X 


CO 


y 


CO 

L z ^ 


(4-25) 


to 


Denoting the Euler angle vector by e_ and the transformation matrix by T(e_) 
which is a function of j:, we get 

* T(e) co (4-26) 

Note that, in (4-25), the transformation equation has singular points at 
8 * mr, n * integer. At these values of 0, sin0 ■ 0 and s ^q does not exist. 

For a strapdown navigation system, the body fixed sensors measure the 

rates w , co and co . The (4-26) is used to generate the Euler angle rates 
x y z 
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Consider a vector V in space. Another vector V' is obtained by 
rotating V about an axis in n direction through an angle 0 as shown in 
Fig. 4-7. Let 1 be the unit vector in n direction, which is related to 
a cartesian coordinate system by 


AAA 


1 “ ^,1+ n 0 j + n. k (4-27) 

n i i -J 

where i, j, k are the three unit vectors of the cartesian coordinates, 

and n, > n-j ru are the three direction cosines of 1 with respect to 

* „ 11 2 3 n 

i, j, k as shown in Fig. 4-8. 


and 



Figure 4-8. Relationship Between i, j, k and 1 . 


From Fig. 4-7, 

V - (1 • V) 1 + |l x V| 1 


n 


(nxv)xn 


(4-28) 


V’ - (1 • V’) 1 + |l x V’ I cos0 1 

n n 1 n 


, . + 1 1 x V’ | sine 1 

(nxv)xn 1 n ' nxv 


(4-29) 


Adding and subtracting (1 q * V) cos6 1^ to (4-29), 


v - u n • v’)(i - cose) i n + (i n x v'| sine i nxv 


+ [(1 • V) 1 + |l x V'l L . ] cos6 

1 n n 1 n 1 (nxv)xn 

N— 


(4-30) 


V 
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But 


|l x V' | 
1 n 1 


1 

nxv 


|l n x V | 


1 

nxv 


1 x V 
n 


1 • V’ •= 1 • V 

n n 

Equation (4-30) becomes 

V* = (1 • V)(l - cos0) 1 + V cos0 + (1 x V) sin0 

n n n 

Adding and subtracting V, 

V’ « V + (1 x V) sin0 + (1 • V)(l - cos0) 1 - (1 - cos©) V 

n n n 

■= V + (1 x V) sin0 + [(1 • V) 1 - V] (1 - cos0). 

n n n 

Using the relationship 

1 x (1 x V) = (1 • V) 1 - (1 • 1 ) V - (1 • V) 1 - V 

nn n nnn n n 

we get 

V* - V + (1 x V) sin© +1 x (1 x V) (1 — cos0) (4-31) 

n n n 

Now let us relate the quantities in (4-31) to the cartesian coordinates and 
express them in the vector-matrix form. 


1 

rH 

> 


< 
M - 
1 


'V 

V_ 

_v ' * 

< 

It 

ho 

2 


2 

n 

2 

. V 3- 


v' 

L v 3 J 


. n 3. 


We then see that 
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”V3 

1 x V = n v 
n J i 

■ n l V 2 

" n 3 V 2 ‘ 
" V3 

“ Vi 

= 

' 0 

n 3 

-" n 2 

" n 3 

0 

n l 

”2 | [ V 

-"l v 2 ■ 

0 J U- 

(4-32) 

where 


0 

“ n 3 

V 




N = 

n 3 

0 

~ n l 


(4-33) 



." n 2 

n l 

0 




is a "direction cosine matrix". Also 

1 x (1 x V) ■= N 2 V (4-34) 

n n 

Hence, (4-31) can be expressed in the form 

V' = EV (4-35) 

where 

E - I + (sin0)N + (1 - cos9)N 2 (4-36) 


is the desired "transformation matrix", which is determined by the angle 
of rotation 0 and the direction cosine matrix N. Since N consists of three 
elements, three are a total of four parameters, namely, 0, r^, and 

They are called "Euler four parameters". Among these four parameters, only 
three are independent, the fourth one can be eliminated by the relationship 


n 


2 

1 




= 1 


(4-3') 


Differential Equations for Euler's Four Parameters 

In the next section, we shall see that the differential equations for 
the quaternion method of coordinate transformation is given by 
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• 



-to 

X 

-to 

y 

-to 

z 



D 


u 

2 

2 

2 


s 

• 

V 

X 


<0 

X 

2 

0 

-to 

z 

to 

y 


V 

X 

• 

V 

y 


tO 

-L 

2 

(0 

z 

0 

-to 

X 


V 

y 

V 

z __ 


to 

z 

2 

-to 

y 

to 

X 

0 


V 

L z J 


(4— 38) 


where s, v , v , and v are the four parameters of quaternion operator. Eq. 
x y z 

(4-38) can also be expressed as 


"* • “ 




s 


-v 

-v 

-v 



X 

y 

z 

• 





V 


s 

V 

-V 

X 

1 


z 

y 

• 

2 




V 


-v 

s 

V 

y 


z 


X 

• 





V 


V 

-v 

s 

w z 


L y 

X 




(4-39) 


It will also be shown in the next section that the Euler's four parameters 
are related to the quaternion parameters by 


0 

cos 2 

. 6 
r \ 1 sin^ 

. e 

n 2 Sltt 2 

, 9 
n 3 sin^- 

Dif ferentiating (4-40), gives 



> 


(4-40) 
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- ~ siny 0 0 0 


•CD 

i . 0 . e A 

2^ sirrj sinj 0 0 


"l 

10 0 
■jr^ 0 sinj 0 


*2 

10 0 
sinj 0 0 sinj 


A 


W 


(4-41) 


Substituting (4-41) and (4-40) to the left and right sides of (4-39), 
respectively, and then premultiplying both sides of the equation by 
yields the desired differential equation 


• 

6 


- 2n l 

-2n 2 


” 2ri 3 ’ 

• 

n x 

l 

Cl - n 1 2 )cot| 

- t, l n 2 cot f + "3 

~ n 3 n l 

9 

cot? “ n 2 

• 

n 2 

“ 2 

• n l n 2 cot f ' n 3 

(1 - n 2 ^)cotj 

-n 2 n 3 

0 

cot^- - n 1 

Ly 


.'"S' 1 ! COC f + "2 

-n 2 n 3 COt 2 " n l 

(1 - 

2, 0 
n 3 )cotj 


>✓ " — > 


u 

x 


0 ) 

y 




(4-42) 


S 


Equations (4-35) and (4-42) form the algorithm for the Euler four 
parameter transformation. The computation requirement for this algorithm 
is shown in Fig. 4-9. 



Figure 4-9. Euler’s 4 Parr er Transformation. 
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QUATERNION ALGORITHM 


Quaternion Algebra 

A "quaternion” Q is defined as 


Q * s + i v + 1 v + k v ■ s + V 
x x J y z 

AAA 

where s is a scalar and V = i v + 1 v +kv is a vector. 

x J y z 

"Addition" for quaternions is defined as 

'l +V 2 


( 4 - 43 ) 


Qj + Qj - (s 2 + s 2 ) + (V, + V,) 


’ (S 1 + s 2> + 1 (v lx + v 2x )+ ^ (v ly + v 2y ) + k (v lx + v 2z> 

( 4 - 44 ) 

It is obvious that addition is associative and commutative. 

"Equality" is defined as 


Ql - Q 2 if only if ■ s 2 & » V 2 


( 4 - 45 ) 


"Multiplication" is defined by using the distributive law on the 
elements as in ordinary algebra, but with the order preserved, and by de- 

AAA 

fining the products of unit vectors i, j, k as 


j 2 \2 :2 , 

i « j = k = - 1 


( 4 - 46 ) 


A A A 


. A A 


ij - - ji = k jk - - kj - i ki - - ik - j 


We can then write 


Q x Q 2 - (S, + v,)(s, + v,) 


1 ' '2 


•l s 2 + S 1 V 2 + S 2 V 1 + V 1 V 2 


V 1 V 2 ‘ (1 v lx + j V ly + k V lz ) (1 v 2x + j v 2y + k V 2z ) 


But 
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- V, V- - V. V- - V. V, 

lx 2x ly 2y lz 2z 


+ i v. v„ - i v, v_ 
ly 2z lz 2y 


+ i V 2x v lz - J v lx v 2z 


+ k v, v_ ■ k v. v- 
lx 2y ly 2x 


- (V, • V 2 ) + (V x X V 2 ) 


therefore 


Q 1 Q 2 " ts 1 s 2 - V x • V 2 ] + [s 1 V^ 2 + s 2 V 1 + V x x V 2 ] (4-47) 


Note (i) the product of quaternions is defined differently from both of the 
dot and the cross products associated with pure vectors, (ii) the product 
of quaternions is not commutative due to the term x in (4-47), and 
(iii) the multiplication is distributive and associative, i.e. 

«1«2 + V * Q 1 Q 2 + Q 1 Q 3 

q 2 (q 2 q 3 ) - <q 3 q 2 ) q 3 

"Conjugate" Q* of Q is defined as 

Q* « 8 _ V (4-48) 

The conjugate of product is the product of conjugate in reverse order, i.e. 

t Q 1 Q 2 ) * - [(s 1 s 2 - V • V) + ( Sl V 2 + s 2 V x + V 1 x V 2 )]* 

- (® X ® 2 - v x • V 2 ) - (s^ + s’ 2 V 1 + V x x V 2 ) 

- QJ Q* 

"Norm" of a quaternion, N(Q), is a scalar defined as 
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N(Q) - s 2 + v x 2 + v y 2 + v z 2 = QQ* 

A unit norm is defined as N(Q) ■ 1. The norm of the product is the 
of the norms, i.e. , 

NCQ^) « (Q 1 Q 2 )(Q 1 Q 2 )* = Q^QJQ* 


Q x N(Q 2 )QJ - N(Q 2 )Q 1 Q* - N(Q 2 ) N(Q 1 ) 


"Inverse" Q ^ of Q is defined as 


Q" 1 Q - qq" 1 - 1 


Since N(Q) = QQ* ■ Q*Q , so 


r 1 .-9 1 


N(Q) 
-1 


For a unit quaternion, N(Q) *= 1, so Q ■ Q* 


Quaternion as Transformation Operator 
Let 

Q “ s + V and N(Q) - 1 

Consider two vectors A and B related by quaternion operator Q via 


A - QBQ* 

Now 

QBQ* * (s + V) B (s - V) 

- (s + V) [B • V + (sB - B x V)] 

- s (B • V) - V • (sB - B X V) + V x (sB - B 

+ s (sB - B x V) + V (B • V) 

- s 2 B + 2sV x B + (B • V) V - V x (B x V) 

2 

Adding and subtracting v B, where 

2 „ „ 2 2 2 
v « V * V ■ v +v +v 
x y z 


(4-49) 

product 


(4-50) 

(4-51) 

(4-52) 

(4-53) 

V) 


QBQ* * (s + v )B + 2sV x B + V x (V x B) 


+ t (B • V)V - (V • V)BJ 

% „ ' 

V x (V x B) 

- (s 2 + v 2 )B + 2sV x B + 2V x (V x B) 

2 2 

from (4-52), s + v - 1, therefore (4-53) becomes 

A - B + 2sV x B + 2V x (V x B) (4-54) 

Expressing (4-54) in the vector-matrix notation in terms of cartesian 
coordinates x, y, z, using the following correspondences 


92 


I 


I 


V x B 

Let 

P - 


(4-57) 

(4-58) 

(4-59) 

(4-60) 


V x (V x B) ♦- 


P 2 b 


Then (4-54) becomes 


Let 


a - [I + 2 a P + 2 P ] b 


W 


W • I + 2 s P + 2 P 


Then (4-58) becomes 


a ■ K b 


0 -v v 

z y 


v 0 -v 

Z X 


, -v V 

L y * 


zj 


(4-55) 


0 -v 


-v 


(4-56) 


-v 


Eq. (4-58) or (4-60) relates vectors a_ and b^ through 
formation W which Is composed of four paramaters s, v 


the matrix trans- 


, v and v . 
x’ y z 
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Isomorphism Between Quaternion and Direction Cosine Transformation 


Quaternion 

Direction Cosine 

U* - Q x U Q* 

u' = C 1 ii 

f = q 2 u’ q* 

u” = C 2 u’ 

>" = q 2 q 1 u q*q* 

u’ ' - C^u 

- (Q 2 Q 1 ) U (Q 2 Q 1 )* 

*= C _u 

*= Q U Q* 


Q - Q 2 Q x 

c = C 2Cl 


Hence, there exists a one-to-one correspondence between direction cosine 
matrices and quaternion transformation elements. 

Equivalence of Quaternion and Euler’s 4-Parameter Method 

Recall (4-53) the quaternion transformation equation. 


a = [I + 2sP + 2P 2 ] b 


(4-61) 


where 


0 -v v 

z y 


v 0 
z 


-V V 

y x 


-V 


0 J 


(4-62) 


Also, recall (4-35), the Euler 4-parameter transformation equation, 


£ * [I + (sin0) N + (1 - cos0) N ] b^ 


(4-63) 
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where 


N 


0 ' n 3 n 2 


"3 ° ‘ n l 


‘-"a 


(4-64) 


Substituting (4-62) into (4-61), we get, for the quaternion method, 


a = < 


t— 1 

o 

o 


r~ 

> 

> 

1 

O 


r 2 2 

-V -v V V V V 



z y 


z y x y x z 

0 10 

+ 2s 

> 

» 

o 

> 

+ 2 

2 2 

V V -v -v V V 



Z X 


x y x z y z 





2 2 

0 0 1 


-v v 0 


V V V V -V -v 


y x J 


x z z y y 3j 


>i 


1 - 2(v + v ) 
z y 

2 (— s v + v v ) 
z x y 

2 (s v + v v ) 
y x z 

2(s v + v v ) 
z x y 

2 2 
1 - 2(v + v ) 
z x 

2(-s v + v v ) 
x z y 

2(-s v + v v ) 
y x z 

2 (s v + v v ) 
x z y 

2 2 
1 - 2(v + v ) 
y x 


(4-65) 


Substituting (4-64) into (4-65) gives for the Euler 4-parameter method. 


a *< 


i o o' 


0 -n 3 n 2 ' 

0 10 

+ sin9 

n 3 0 _r| l 

.0 0 1 


-n 2 n 3 0 


+ (1 - COS0) 


" n 3 ~ n 2 


n l n 2 


n l n 2 n l n 3 


2 2 

‘ n l “ n 3 n 2 n 3 


u n 1 n 3 


n 3 n 2 


2 2 
-n 2 “ \ j 


> ° 
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0 0 ** 0 

Since sin0 = 2siny cosy and 1 - cos6 = 2sin - this equation can be 
expressed as 



1 

0 

0 


0 

, 0 
'V in 2 

, 0] 
V ln 2 

< 

0 

1 

0 

0 

+ 2cosy 

. 0 
n 3 sin 2 

0 

. e 

-n^siny 


. 0 

0 

1 . 


, 6 
-V in 2 

. 0 
n^sin^- 

0 . 


s 


' 4 * 

-sin y 

2 2 
(tl^+np 

. 0 

sin y 0^2 

sin y 

010 3 


1 20 
sin ^ 

V2 

. 2 0 / 2, 2 N 
-sin y (0^03) 

. 2 0 
sin y 

0 2 0 3 


2 e 

l sin 2 

n l n 3 

sinf n 2 n 3 

2 0 

-sm y 

2 2 
(01+ n 2 )J 



Also recall that for the quaternion method, 


(4-66) 


A - B + 2sV x B + 2V x (V x B) 


(4-67) 


and for Euler’s 4-parameter method, 

A = B + (1 x B) sin6 + 1 x (1 x B) (1 - cosO) (4-68) 

n n n 

Let us make the following identification: 

If V is f-’kon as a vector along the 1 direction, then 

n 

2sv = sin0 (4-69) 

2v 2 = (1 - cos0) 

Solving for v and s from (4-69), 

v * t si try j 

+ 0 
8*3 COSy 


> 


(4-7C) 
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Since (4-70) satisfies N(Q) *- 1 irrespective of signs chosen and since 

A *= QBQ* = (-Q) B (-Q)* (4-71) 

the rotation quaternion is arbitrarily chosen using positive signs in (4-70) 

Q ■ s + V = cos^- + 1 sirr^ ■» s + 1 v (4-72) 

2 . n 2 . n 

where 


6 

s « cos^- 


v » sin^ 


(4-73) 


Comparing (4-64) and (4-66) and using (4-73), we have the following equivalence 
between the quaternion and Euler’s 4-parameter methods: 


s = cos^- 


V = 
X 

cosa 

. 0 

sin^- - 


. e 
sln 2 

< 

II 

COS0 

. 0 

Sin 2 = 

n 2 

. 6 
sirr 2 

V = 

z 

cos y 
\ 

. 0 

sm-j = 
y"_ 

"3 

, e 

sin^- 


(4-74) 


Quaternion 

Method 


Euler’s 4-Parameter 
Method 


Parameters: 
s » v „> v ..» v . 


Parameters: 

x y z 12 3 

The relationship among various variables is shown in Fig. 4-10. 
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Figure 4-10. Relationship Between a, B, y and 1^. 

Among 4 quaternion parameters, only 3 are independent; the 4th one can be 
eliminated by using the constraint N(Q) = 1. 


laternion Differential Equations 


Let Q(t) be a quaternion representing the transformation from a rotating 
body frame to the inertial frame. Let be the value of Q(t) at time t, i.e.. 


A i (t) = QjA b (t) Qf (4-75) 

and Q 2 be the value of Q(t) at time t+At, i.e., 

A 1 (t+At) = Q 2 A b (t+At) Q* (4-76) 

The quaternion representing the rotation of the body frame during the time 
At is therefore, 

qTV - Q*Q 

x z x 2 

This relationship is shown in Fig. 4-11. 



Figure 4-11. Quaternion snd Vector Rotation. 
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Since the angular velocity of the rotating body frame can be considered 
constant for small At the angular displacement of the body frame is 

A0 - loAt, 


where 


n 


- and | d | 


Therefore 


and 


qjq 2 


COJ 


wAt 


ft , wAt 
+ — sin-r- 
u) 2 


(4-77) 


_ . F QjAt ft . toAt 

Q 2 " Q 1 t cos — + u sln — 


(4-7 b) 


The time derivation of the quaternion Q is obtained from 
Q (t) 


lim Q 2~ Q 1 

lim Q_ 

At+o At 

At->o At 

lim Q f 
At->-o At j. 

i - > 2 

1 


f Qfi 


1 < 8 + v > 



coAt , , ft 0)At 

cos— r— - 1 + — 

2 0 ) 


- , ft /OiAtv 1 

- 1 + * - — J 


• ft) + (s ft + 


V x ft) ] 


(4-79) 


Hence, we get the desired differential equations 


8 - - Y [V • fl] 


v- j [s ft + v x ft] 


To express in vector-matrix rotation, let 


s 

v 


Q * (s,V) - 


(4-80) 

(4-81) 


(4-82) 


z J 


99 


Then (4-80) and (4-81) can be expressed as 


■ . « 
s 


• 


V 


X 


• 


V 


y 


• 


N 

> 

t 




U) 

0) 

0) 

o 

- X 


- z 

V/ 

2 

2 

2 

0) 


0) 

0) 

X 

o 

- z 


2 


2 

2 

w 

0) 


U) 

JL 

z 

0 

- X 

2 

2 


2 

0) 

0) 

0) 


z 


X 

o 

0 


(4-83) 


v 

*• z 


M 


A 1 (t+At) = Q 2 Q b A b (t) Q b Q* 

- Q 2 Q*Q* A*(t) Q 1 Q b Q* 


Since V*(t+At) = V*(t) , the stationary inertial frame, 


Wi = 1 

Q b - QJQ 2 (4-84) 

The computation procedure from body axis measurements to inertial axis 
measurements is shown in Fig. 4-12. 



Figure 4-12. Quaternion Transformation. 
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Remarks 

1. Transformation of coordinates can more generally be expressed as 
quaternion operated on quaternion. Let v and ju be the same 3-dimensional 
vector represented WRT two different coordinate systems. Let 


Q 8 + V = operator quaternion 

■ s^ + V = quaternion representing vector v 

Q 2 * 82 + ■ quaternion representing vector u 

Then (4-53) can be expressed as 

Q 2 - QQ 1 Q* 

- (s + V) (Sj, + VjKs - v) 

- (s + V) s L (s - V) + (s + V) V x (s - V) 

1 ■ ■"’* * 1 ■ ■ 1 111,1 ^ 

8 ^ let it be V 2 


«° s. + V 

X 


2 


(4-85) 


Thus, the operation does not change the value of the scalar term. 

2. Let represent the rotation of the body frame from time t to t+At 


A b (t+At) « Q* A b (t) Q b (4-86) 

Then at time t 

A*(t) = Q A b (t) Q* (4-87) 


and at time t+At 

A 1 (t+At) 
Since A 1 (t+At) = A i (t) 


A 1 (t+At) - Q 2 A b (t+At) Q* 

- Q 2 Q* A b (t) Q b Q* - Q 2 Q{Q 1 A*(t) Q 1 Q b Q 2 
, the stationary frame vector, so Q 2 Q£Qi 


(4-88) 


1 . Hence 


Q b - QJQ 2 
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CAYLEY-KLEIN PARAMETER ALGORITHM 


Transformation In Two Dimensional Complex Space 

Consider the transformation of two complex axes (u,v) to (u',v') 
through the linear relationship 



(4-89) 


Where a, 8, y and 6, called Cayley-Klein Parameters, are all complex numbers. 
Eq. (4-89) denotes a transformation in a two-dimensional complex space, or, 
equivalently, in a four-dimensional real space. 

Eg. Transformation in 1-dimensional comp.i ;x 
space 

(a’+jb*) •= (a+j B) (a+jb) 

can also be expressed as transformation 
in 2-dimensional real space as 

[::] ■ [: -.] [:] 

Since a, B, y and <5 are complex, there are eight parameters to be specified. 
However, if one tries to apply the K matrix to a transformation in a 3- 
dimensional real space, only three of the eight parameters are independent. 

To eliminate the remaining five parameters, some constraints are needed. 

The constraints are: 

1. K be a unitary matrix, i.e., K * ■ K*, or just 

K*K « 1 (4-90) 

2. The determinant of K be 1, i.e., 


K| - 1 


(4-91) 
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Using (4-89) 


requiring 


From (4-91) 


K*K 


a 

Y 


a* 

Y* 

8* 

6* 


a*a 

+ 8*8 

aY* 

+ 66* 


» 

H* 

O 

1 

a*y 

+ 8*6 

YY* 

+ 66*_ 


-° 


oa* + 66 * 

■ 1 

(4-92) 

y-f + 6 6 * 

- 1 

(4-93) 

a* y + 6*6 

- 0 

(4-94) 

06 - 6 Y 

- 1 

(4-95) 


Eqs. (4-92) and (4-93) are real, while (4-94) is complex, so that together 
they provide four conditions. Eq. (4-95) provides the fifth condition. 
From (4-94) 


& 

Y 


a* 

6 * 


(4-96) 


Substituting into (4-95) 


ua* 



36* (-^r) ■ 1 


(aa* + 88*) * - 1 


By the relationship of (4-92), we have 

Y - - 6 * (4-97) 


Using (4-97) in (4-96) 


6 - - a* 


(4-98) 
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Therefore, the matrix K can be written as 


Constrained by 


a 


K - 



8 



aa* + 88* * 1 


(4-99) 


(4-100) 


Notice that the derivation of (4-99) used (4-92), (4-94), (4-95) but not 
(4-93). (4-93) is used as the constraint to (4-99). On account of (4-97) 

and (4-98), (4-93) ■ (4-92), so we have (4-100). 


Now let D be a Hermitian matrix in this 2-dimensional complex space 
having the specific forr 


D - 


x 

z+jy 


z-jy 

-x 


J - 


(4-101) 


The quantities x, y, z are reals, and they can be considered as representing 
three coordinates defining a vector r in a three dimensional real space. 
Notice that the determinant of D is 


| D [ - - (x 2 + y 2 + z 2 ) = -r T r «= - r 2 (4-102) 

and 

Trace D - 0 (4-103) 

Applying a similarity transformation to V using the unitary matrix K, gives 

D' - K D K* (4-104) 

Under similarity transformation the Hermitian property of a matrix is pre- 
served, so are the determinant and trace, i.e., 

|D* | “ “ (x’ 2 + y ,2 + z' 2 ) - | D f 


Trace D' - 0 
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Isomorphism Between Cayley-Klein Parameters and Quaternion Parameters 
Matrix D can be written as 


D - 


x 


L z+Jy 


j-jy 



‘l 0 ' 


'o -j ' 


‘o 1* 

■ X 


+ y 


+ z 



0 -1 
• ■ 


O 

*n 


.1 o. 


■ x E + y I + z 2 
x y z 


where 


In addition 



(4-107) 


(4-108) 


Since these 4 matrices are linearly independent, any 2x2 matrix involving 
only 4 Independent quantities can be formed by linear combination of them. 
Let 

a * a + jb 


8 ■ c + jd 


> 


(4-109) 
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Then using (4-108) we have 
' a 8 

->8* a* 


a+jb 

c+jd 

-c+jb 

a-jb 


a< 


"l o‘ 


. 

1 0 


'o -j ' 


' 0 1 ' 


+ j b 


+ c 


+ d 


.0 1. 


.0 -1 . 


J 0. 


t-* 

o 


A K = a E T + j I bE + cZ + dZ I 
I J l x y zj 

Recalling the quaternion Q, 


Q *» s + V •= s + iv + 1v + kv 
^ x J y z 


(4-110) 


(4-111) 


There is an isomorphism relationship between K and Q where: 


K ► Q 

a *- s 



d 

V 

z 


j£ x i 

i 




J 


3£ x ► 

A 

k 

Addition: 

h + *2 - 

Ql + Q 2 

Equality: 

K 1 * K 2 

- **2 


Multiplication: 
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a y 2 - ay 2 - ay 2 - - 


A rt A A 

J - k" 


(JS x )(jS y ) - - <jz y )(j£ x > - 


- ji = k 


<3 E y )( JE «) = ' <jE z )(J V = jE x 


ay ay - - ay a= z ) = 


- kj - i 


- Ik = j 


K 1 K 2 " ta l E I + b l (jE x ) + c l (jE y ) + d l (j V ] 


[• 2 Sj + b 2 ay + c 2 ay + d 2 ay ] 


[ V2 " (b l b 2 + C 1 C 2 + d l d 2 )] E I 


+ * x [b 2 ay + c 2 ay + yjyi 
+ a 2 [b x ay + c 1 ay + yjy i 


IE jE jE 

J x J y J z 

+ bj c i d i 


b 2 c 2 d 2 


Q 1 Q 2 ’ <S 1 + V (S 2 + V 


(8 1 + iv lx + Jv ly + kv l* )(s 2 + lv 2x + JV 2y + kv 2x ) 


(s l s 2 ' V lV + (S 1 V 2 + S 2 V 1 + V 1 x V 


b l b 2 + c l c 2 + d l d 2 


V • V 
1 2 
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»lIb 2 OS x ) + + d 2 < J E « )] S 1 V 2 

a 2 [b l Cj£ X > + c l^^y) + d 1 (jE z )] „ ». SjVj 

JI X J£ z 

b l C 1 ^ 

b 2 c 2 d 2 



| Norm: |k| *= a^+ b^+ c^+ *■ N(Q) = s^+ v^+ v^+ 

| That is, in Cayley-Klein parameters, the norm is taken to the determinant 

of K matrix. 

Conjugation: K* = transpose conjugate of K « -**Q* = s - V 

Transformation of Coordinates 

In view of the isomorphism between the Cayley-Klein algebra and the 
quaternion algebra, and in view of the fact that vector transformation in 
quaternion algebra is represented by 

Q 2 - QQjQ* (4-112) 

therefore, the vector transformation in Cayley-Klein algebra can be expressed 
as 

D 2 - KDjK* (4-113) 

where the two sets of parameters, in the operators K and Q, are made identical 
to each other in the corresponding terms. That is, 

Cayley-Klein Quaternion 

a « s 

b ■ . v x 

I C - Vy 

d “ v z 

I 
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Therefore, basing on the result of quaternion transformation method, we can 
express the Cayley-Klein transformation method as 


where 


u 



[1 + 2aP + 2P 2 ) v 
G 

-d c 

0 -b 

b 0 . 


(4-115) 


(4-116) 


Cayley-Klein Parameter Differential Equations 


Due to the relationship (4-114) we have, based on the differential 


equations for quaternions, 



(4-117) 


M 


The computation procedure from body axis information to inertial axis 
quantities is shown in Fig. 4-13. 



v 


Figure 4-13. Cayley-Klein Parameter Transformation. 
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VECTOR REPRESENTATION ALGORITHM 
Vector Representation 

Let 1 be the unit vector along the axis of rotation, and 0 be 
n 

the angular displacement as shown in Fig. 4-14. 



Figure 4-14. Vector Representation. 


Define a vector E 


1 tanj- 
n 4 


. 0 
sin-7 

^n-r 

cos 4 


n 


. 6 0 

sin-r cost 
4 4 

2 0 
COS *7 


, 0 
siltj 


. 0 
sinj 


o 0 

2 cos ~jr 


n . , 0 

1 + cos^- 


Recall that in the quaternion method where Q = s + V, we have 

6 6 
s *= coar- and V = 1 sinrr 

2 n L 


Therefore 


1+s * 


(4 


-119) 


so 


V - (1 + s) E 


(4-120) 
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Solve for s and V in terms of E . Note 

2 


- , . . 9 _ 9 2 

1 + s = 1 + COST- = 2 cos -r •= r 

2 4 2 


0 

CSC T 

4 


1+tan y 
4 


9 1+0 


2 2 

where o *» £ ■ tan -r 

4 


1-g 

1+0 5 


(4-121) 


Substituting (4-120) into (4-119) 


l+o 


Recalling the quaternion transformation equation, 


A - B + 2s (V x B) + 2 V x (V x B) 


(4-122) 


then by using (4-121) and (4-122), we have the transformation in using the 
vector representation method. 


A = B + 4^‘V - , (E x B) + 

(l+o 2 ) 2 (l+o 2 ) 2 


E x (E x B) (4-123) 


Expressing in vector-matrix rotation 


E x B 


-o 


-o 

«- y 


-o 


. 


lb 2 j 


(4-124) 


and 


E x (E x B) « S*b 


(4-125) 


where 


A ^ ^ 

B-ib + j b + k b 

x J y z 





Ill 


Using (4-124) and (4-125) in (4-123) 


-[ 


! + 4 U-O s t - 8 

(i +o (i .+o z 


* 


T 


(4-126) 


Let 


T « I + 4 - ^~ a } -■ S + — 

22 22 

(l+o ) (1+cr ) 


(4-127) 


(4-126) becomes 


a ■ T b 


(4-128) 


Differentiation Equation for Vector Representation 
Differentiating (4-120) 


• • 


V » (1 + s) E + s I 


(4-129) 


Since 


s - — f V • ft 


V- ~ (s 0 + V x 0) 


(4-129) gives 


(1 + s) Z - j (s Q + V x 0) + I (V • fl)Z (4-130) 


Using (4-120) in (4-130) 


(1 + s)E - j [s 0 + (1 + s) (E x fl) + (1 + s)(Z • 0) S] 

^“l[l+s fi + (E x 0) + (E • fi) e] (4-131) 


From (4-121) 


s 

1+s 


1-0 
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Using this in (4-131) 


1 f hsl n 

2 [ 2 


+ (r x ft) + (E • n)E 




] 


Nonlinear 
terms 


By using the identity 


we have 


A x (B x C) - (A • C)B - (A • B)C 


r x (E x ft) - (i • n)E - (2 • z)n 


(E • ft)E - oft + Z x (Z x ft) 


(4-132) 


Therefore (4-132) can also be written as 


f 1 {* l+o 2 

[ZJZ 


SI + (L x S!) + l x (l x S!) (4-133) 


] 


V 


Nonlinear 

terms 


DDA REALIZATIONS OF ATTITUDE ALGORITHMS 

Direction Cosine Algorithm 
The transformation equations are 



C-Cl 



(4-17) 


2 1 - C V B 



(4-18) 

Taking the differentials and 

let ait - 4A, 

the angular increments. 


4C u 

■ C i2**3 - 

C 13 “2 ' 



4C i2 

' °13 4A 1 - 

C 11 4A 3 

► i - 1,2,3 

(4-134) 

4C 13 

- C 11 4A 2 - 

C 12 M 1 . 
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AV 


AV 


AV 


C 11 C 12 


C 21 C 22 


C 31 C 32 


'13 


'23 


AV 


AV 


'33 J 


AV 


3 J 


(4-135) 


The DDA realization for (4-134) is shown in Fig. 4- 15a. The realization 
for all i - 1,2,3 requires 6 x 3 ■ 18 DDA’s. The DDA realization for (4-135) 
is shown in Fig. 4- 15b. The realization for all i « 1,2,3 requires 3x3“ 

9 DDA’s. Therefore, the complete realization of both (4-134) and (4-135) 
for all coordinate variables requires a total of 18 + 9 ■ 27 DDA’s. 


Euler Angle Algorithm 

The Euler angle differential equations are given by 


In incremental form. 


A<p ' 


A0 

- 

Aip 

► . i 



sinip 

cosip 

o 

sin6 

sin0 


cosip 

-sinip 

0 

sinip cos0 

COS0 cosip 


sin0 

8in0 

1 

si nip 

CQ81p 

o 

8in0 

sin0 

V 

cosip 

-sinip 

0 

sinip cos0 

cosip COS0 

1 

sin0 

8in0 

i. 


u> 


a) At 
x 


a) At 

y 


0) At 
z 


(4-25) 


Let 


AA X - uj^At 


AA^ “ At 


AA- « oj At 
j z 


(4-136) 
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We get 


- sini/> + AA 2 cos)iJ/ — j^- 


A0 « AA^ cost p - A A 2 sin t|> 


A>J> - (AA^ sinif* + AA 2 cos \p ) — —q cos9 + AA^ 

The transformation equations are given by 

cos4< cos<f» - cos0 sin$ sin^ 

-sin^ cos$ - cos0 sin<J> cosi p 
sin0 sin<f) 


< 

' 


m 

.< 



cosi^ sin<j> + cos0 cos<j> sinij j 
- eitvfi sin4> + cos0 cos<J> cos<J» 
-sinP cos<f> 


sirup sin9 


'< 

C 08 ip sin0 



cos0 

m 


< 


Furthermore, the following relationships exist. 


Asin$ « 

A$ 

cos^i 

Acos$ * - A$ 

sin$ 

Asln0 * 

A0 

cos$ 

Acos6 ■ - A0 

sin0 

AsiniJ> ■ 

Aip 

cos^ 

Acost|> “ - A^i 

siniji 



A — - 

AsinO 




sin0 

sin Z 0 



(4-137) 


(4-138) 


(4-139) 


The DDA realization of (4-137) is shown in Fig. 4-16a» using 6 DDA's and 
the realization of (4-139) is shown in Fig. 4-16b> using 8 DDA's. The DDA 
realization of (4-138) is shown in Fig. 4-16ct using 13 DDA's. Therefore, 
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the complete realization of all (4-137), (4-138) and (4-139) for all 
coordinate variables requires a total of 6 + 8 + 13 = 27 DDA's. 

Euler’s 4 Parameter Algorithm 
The Differential equation is 


’ e ‘ 


r - 2 n 

-2n 

-2n 


0) 



X 

y 

z 


X 

• 



e 

0 



n 


cot-r- (1-r ) 

~n n cotr- -n 

-n n cot r- +n 



X 

l 

2 x 

x y 2 z 

x z 2 y 




" 2 





(1) 

• 


e 


0 


y 

n 


-n n cotr- 

cotr (1-n ) 

-n n cot-r- -n 



y 


x y 2 z 

2 y 

z y 2 x 



• 


0 

0 

0 2 n 



n 


-n n cot x- -n 

-n n cotrr +n 

cot-r- (1-n ) 


0) 

z 


X z 2 y 

L 

y z 2 x 

2 z J 


z . 


or, in the incremental form. 


ri 


-2n 

X 

-2n 

y 

-2n 

z 

> 

rt 

X 

. * 

.0 2, 
coty Cl-n x ) 

0 

-n n cot-r- -n 
x y 2 z 

0 , 

-n n cotr- +n 
X y 2 y 

A "y 

JL 

0 , 

-n n cotr- +n 
X y 2 y 

cot^- (l-n ) 

0 

-n n cot-r- -n 
z y 2 x 

, 4n z. 


0 

-n n cotrr- -n 
L x z 2 y 

e 

-n n cot^r +n 

y z 2 x 

6 2, 
coe 2 (1_n z ) J 


AA, 


AA, 


AA, 


(4 


The transformation equation is 


AV 1 = [I + N sin0 + N 2 (l - cos0)] AV® 

where 



-140) 
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Or, we can write it in the form of 


AV 1 *= AV B + sin9 1 x AV B + (1 - cos6) 1 x (1 x AV B ) (4-141) 

n n n 


1 »in + j n tk n 
n x J y z 


Eq. (4-141) can be written in the following more explicit form 


Av^ “ Av B + sin0 • + (1 


X B 

AV 2 •= AV 2 + sin0 • + (1 


Av* = Av B + sin0 • U 2 + (1 


cos0)(u„ti - u.n ) 
Z X I z 


cos0)(u 3 n z - u 2 n x ) 


cose)( U;L n x - u 3 n y ) 


(4-142) 


where 


U 1 - ", 4v l - n * 4v 3 


u. 


"x 4v 2 ' "y 4v l 


"y Av 3 ' ", 4v 2 


(4-143) 


It takes 2 DDA’s to take care of the factor coty using the relationship 


2 2 
Acoty * - CSC y • A0 = - (1 + cot y) A0 


(4-144) 


The realization of (4-144) is shown in Fig. 4-17a. Six DDA's are needed to 
realize (4-143) as shown in Fig. 4-17b, and fifteen are needed to realize 
(4-140) as shown in Fig. 4-17c. It takes 2 DDA's to get the factor sin0 and 
cos0 using the relationships 


} 


(4-145) 


A sin0 = cosG • A0 

-A cos0 ■ sin0 • A0 
The realization is shown in Fig. 4-l7d. Twelve DDA's are needed to realize 
(4-142) as shown in Fig. 4-17e. Therefore, it takes a total of 2 + 15 + 6 + 
12 + 2 - 37 DDA’s. 


CD | CM 
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ler 4 Parameter Algorithm. 
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Quaternion or Cayley-Kleln Algorithms 

The DDA requirement for quaternion method or for Cayley-Klein method 
is the same due to the isomorphism of these two methods. The quaternion 
differential equations are expressed by 


co 


CD 


—CO 


CD 


-CD 


-CD 


-CD 


0 ) 


-CD 


-CO 


(4-83) 


In incremental form 
As 
Av 


Av 


Av 


0 

“i 
4 A„ 


-AA, -AA. -AA, 


0 

AA, 


AA_ -AA_ 


-AA, 

0 

AA, 


AA- 


-AA. 


(4-146) 


where 


AA, 


AA- 


The transformation is given by 

V 2 

cl 


CD At 
X 


CD At 

y 


AA- ■ cd At 
3 Z 


2 2 
1 - 2(v + v ) 
z y 

2(-sv z + vv y ) 

2(sv + v v ) 
y x y 

2(sv + v v ) 
z x y 

1 - 2(v 2 + v 2 ) 
z X 

2(-sv + v v 

x z y 

2(-sv + v v ) 

y x z 

2 (sv + v v ) 
x z y 

2 2 

1- 2(v + v ) 
y x 





V 1 


B 


V 2 


B 

• 

JJ 
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Multiplying the right hand side out and rearranging terms, we get 


Let 




» . * i / jJ Bv _l/B B * , * B B N 

V. + 2 (v 0 v - v,v )v„ + (v.v - v_v )v + (v,v - v 0 v )s 


3 x 1 z z 


2 * v-v/v » \ r « v ▼ a » j 

x 1 y y 3 y 2 z 


] 


v 2 ■ v 2 + 2 


v 3 " v 3 + 2 


[ (v?v - v?v )v + (v?v - v®v )v + (v?v - V^v )s ] 
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Then, we have, from (4-147) 
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The complete coordinate transformation is given by (4-146) , (4-148) and 
(4-149). The DDA realization of (4-146) requires 12 DDA's as shown in 
Fig. 4-18a, the realization of (4-148) needs 6 DDA's as shown in Fig. 

4- 18b, and the realization of (4-149) requires 9 DDA's as shown in Fig. 
4-18c* Therefore, the complete realization of (4-146), (4-148) and (4-149) 
for all coordinate variables requires a total of 12 + 6 + 9 ■ 27 DDA's 
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Figure 4-18. (Continued). 
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Vector Representation Method 
Recall the equation 
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Expanding (4-133) and then to put it into the matrix form 
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Then (4-150) becomes 
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The transformation equation is given by 


V 1 = V B + 4(1 ~° l (Z x V B ) + 8 V 2 (I x Z x V B ) 

(i +a z r (1+0 


(4-153) 


Recognizing the similarity between (4-153 and (4-133) in their form, we can 
write (4-153) directly into the following form using incremental. 
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I B 
Av„ = Av„ 



A (1+0 ) = A (1+0 + 0 + 0 ) = 2(0 Ao +o Ao +0 A 0 ) 

xyz xxyyzz 

it needs 3 DDA’s to realize as shown in Fig. 4-19a. The quantity 


( 1+0 


(l+o ) 


(1+0 ) (l+o ) 


can be realized by 2 DDA's as shown in Fig. 4- 19b. Equations (4-151), 
(4-152), (4-154) and (4-155) can be realized by 6, 9, 12 and 6 DDA's, 
respectively, as shown in Figures 4- 19c, d, e and f. Therefore, a total 
of 3+2+6+9+12 +6= 38 DDA's are required. Notice that 

A(l-o^) = - 2 (0 Ao + 0 Ao + 0 Ao ) = - A (l+o^) 
xxyyzz 

This operation is also included in Fig. 4-19. 


Summary 

The following is a table summarizing and comparing the numbers of 
required DDA's for each attitude algorithm, as obtained from this study 
and obtained by other research.^ - 


Table 4-1. Summary of DDA Requirements. 


Algorithm 


Number of DDA's Required 


Reference 1 This Study 


Direction Cosine Algorithm 
Euler Anple Algorithm 
Euler 4-Parameter Algorithm jj 

Quaternion or Cayley-Klein Algorithms 
Vector Representation Algorithm I 
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CHAPTER 5 


REBALANCE ELECTRONICS 


T.V. Blalock and E.J. Kennedy 


ABSTRACT 


Two basic types of strapdown gyroscope rebal- 
ance-electronics were analyzed and compared. These 
two types were a discrete-pulse ternary system devel- 
oped at the MIT Draper Laboratory and a width-modulated 
binary system developed at the Hamilton Standard System 
Center. In the analyses, major emphasis was placed on 
the logic sections, the H-switches, the precision vol- 
tage reference loops, the noise performance, common- 
mode rejection, and loop compensation. Results of the 
analyses were used in identifying specific advantages 
and disadvantages of system details and in making ac- 
curacy and resolution comparisons. Sound engineering 
principles were applied in the development of both sys- 
tems; however, it was concluded that each system has 
some disadvantages that are amenable to improvement. 


INTRODUCTION 

A preliminary study of three different realizations of rebalance 
electronic loops for strapdown gyroscopes resulted in the selection of two 
of the three for the detailed analyses presented in this chapter. The 

three systems were (1) the Honeywell ternary^, (2) the MIT Draper Labora- 

2 3 

tory ternary , and (3) the Hamilton Standard width-modulated binary . 

After the preliminary study, the Honeywell system was eliminated 
as a candidate for the detailed analyses for two reasons, (1) the 
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preliminary study results indicated that the Honeywell system had several 
disadvantages not apparent in the other two systems, and (2) the allotted 
contract time would allow a detailed analysis of only one example of the 
two basic types of systems — ternary and binary. 

In the second section of this chapter a description of the two 
systems analyzed will be given. In the third section, analyses of the 
logic blocks, H-switches, PVR loops, noise performance, common-mode re- 
jection and compensation are presented. Resolution and accuracy are con- 
sidered in the fourth section. The last section contains a summary of 
comparisons between the two systems analyzed and concluding recommendations. 


DESCRIPTION OF SYSTEMS ANALYZED 

A general description of the organization and operation of the 
Hamilton Standard and the MIT rebalance electronics will be presented in 
this Section. 

Hamilton Standard System 

The basic components of the rebalance electronics for the Hamilton 
Standard width-modulated-binary strapdown gyroscope system is shown in Fig. 
5-1. The system can be divided into two sections, (1) the rebalance loop 
that surrounds the gyro, and (2) the digital control electronics (DCE) which 
generates the basic commands for system operation. 

The rebalance loop senses movement of the gyro float from its quies- 
cent position and then generates pulses of torque current to restore the 
float to its quiescent position. The DCE determines the error signal sam- 
pling frequency, determines the torquing mode and switches scale factor 
when necessary, and generates 256 kHz data pulses which give the amount of 
torque-current area (ampere-seconds) fed into the torquer. 

The input to the rebalance loop is a 50 kHz sinusoidal error signal 
from a sensitive (1. 39mv. /arc-sec) capacitance bridge pick-off mounted in 
the gyroscope. The error signal is processed by a high-input, impedance 
preamplifier which removes common-mode signals contaminating the error 
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signal. The output of the preamplifier feeds a bandpass filter which 
further improves the signal- to-noise ratio of the amplified error signal. 
The bandpass filter feeds a synchronous demodulator whose dc output vol- 
tage can have positive or negative polarity depending on which of the two 
possible directions the gyro axis has moved. The dc error signal is am- 
plified, compensated for system dynamics and then mixed with a 1 kHz ramp 
signal. The zero-voltage crossing point of this composite signal deter- 
mines how the total torque applied to the torquer will be apportioned be- 
tween positive and negative torques. The zero crossing is detected by a 
group of logic gates in the quantizer (Fig. 5-1). Detection of the zero 
crossing is indicated by a transition in the quantizer J-K flip-flop which 
feeds torque command signals to the H-switch drivers and a 1 kHz forced 
limit-cycle signal to the DCE for torque mode determination. The frequency 
(1 kHz) of the limit cycle is set by the ramp reset signal which is the 
inversion of the 241-256 blanking signal generated in the DCE. The blank- 
ing signals (0-16 blank and 241-256 blank) assure that for any torquing 
condition, both positive and negative torque will be applied for at least 
16/256 msec. 

The H-switch routes a precisely controlled current I to the torquer. 
The bilateral nature of the H-switch allows current to be driven in both 
directions through the torquer so that the net movement of the gyro float 
axis for each 1 msec torquing interval is proportional to the difference 
between the positive torque current area and the negative torque current 
area. Thus, a zero net torque requirement is fulfilled by equal positive 
and negative torque current areas. 

The system scale factor is determined by the current level feeding 
the H-switch; consequently, a scale factor change is executed by switching 
the current level with a scale— factor change command from the DCE to the 
torque current regulator. The current level is changed in the regulator 
by switching the values of the current sensing resistors. The current I 
is about 143 ma in the high mode and about 71.5 ma in the low mode. The 
scale factor change command also changes the effective gain in the re- 
balance loop by varying the ramp signal slope. 
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The digital control electronics is fed by a A. 096 MHz clock signal. 

The DCE consists of two primary blocks; the digital function generator (DFG) 
and the dual scale-factor select logic (DSFS). The DFG generates (1) a 256 
kHz, 1 ysec wide clock signal for logic synchronization and data output, 

(2) a 0-16 blanking signal, (3) a 241-256 blanking signal. (A) a ramp reset 
signal, and (5) a set of six mode-interrogation signals. All of these DFG 
signals are generated in each 1 mstc interrogation period. Two of the mode 
interrogation signals determine the allowable high-mode limit-cycle range, 
two determine the low-mode range, and two are used to set the DSFS for scale 
factor change. 

In the DSFS, the mode Interrogation signals examine the 1 kHz limit 
cycle from the quantizer to ascertain in which .torque mode (high or low) the 
rebalance loop should be operating. If the high mode is needed, the DSFS 
will generate a high scale factor command after A high torque-current pulses; 
if the low mode is required, the low scale factor command will be generated 
efter 20A8 low torque current pulses. 

The DSFS also generates the 256 kHz data pulses which contain the 
torquing information. The DSFS selects one data line for high mode, and 
another for lov mode. High and low mode sync pulses are transmitted out on 
two other lines to indicate the initiation of a data pulse train (a data 
train is produced for each 1 msec interrogation period). 


MIT Draper Laboratory System 


The MIT rebalance electronics is a ternary system that is composed of 
a rebalance loop and a logic block (Fig. 5-2). Pulses of torque current are 
routed to the torquer only when the error signal exceeds a pre-set positive 
or negative limit. The polarity of the error signal depends on which of the 
two possible motions of the gyro float axis has occurred, and therefore de- 
termines the polarity of the rebalance torque to be applied. The error sig- 
nal is sampled at 104 usee intervals (9.6 kHz rate). In each 104 ysec in- 
terval, a 6.5 usee period is used to set up the torquing circuitry, and a 
97.5 usee torque command pulse is generated. The torquing circuitry is set 
up to (1) route the torque pulse to the torquer or to a dummy load, and (2) 
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to generate a negative or positive torque pulse. 

The input of the rebalance loop is a 9.6 kHz error signal from the 
gyro float position pickoff which is fed by a 9.6 kHz signal generated in 
the logic block. Thus, the error signal phase is synchronized with the 
command signals generated in the logic block. 

The error signal is first processed by the amplifier and compensa- 
tion section (Fig. 5-2). In this section the error signal is amplified, 
filtered to reduce common-mode and electronic noise, and mixed with a paral- 
lel compensation signal. The compensation is applied only when the torque 
command pulse is present and causes that part of the rebalance loop from 
torquer input to amplifier output to have a response that closely approxi- 
mates that of a pure integrator. The scale-factor command changes the gain 
of the amplifier by switching the value of a feedback resistor in an oper- 
ational amplifier. 

The amplified and compensated error signal (a 9.6 kHz sinusoid) is 
next fed to the comparator. At the comparator input the error signal is a 
sinusoid superimposed on a 2.5 volt dc level. The comparator reference 
signal (see Fig. 5-2) is also 2.5 volts. Since each interrogation period 
spans one cycle of the 9.6 kHz error signal, the polarity of the peak am- 
plitude of the error signal at a fixed strobe time, which is set by logic 
synchronizing signals, determines the polarity of the required torque cur- 
rent pulse, and the amplitude of the error signal determines the need for 
a torque current pulse. That is, if the peak amplitude of the error signal 
is outside the comparator set points, 2.5 + 0.3 volts, a torque current 
pulse of appropriate polarity will be routed to the torquer . The two out- 
puts (one for positive and one for negative error signal polarity) are in- 
terrogated in the logic block by the strobe signal. 

Connected in parallel with the input to the comparator is a sample- 
and-compare circuit. A sample command pulse from the logic block opens a 
field-effect transistor (FET) gate at the time of occurrence of the error 
signal peak. The peak value is then stored on a capacitor at one input 
to a comparator. The other input to the comparator is fed by a staircase 
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signal from the logic block. The number of steps necessary to reverse 
the comparator output is recorded in the logic block which generates a 
number of data pulses proportional to the number of stair steps, and thus, 
performs an analog-to-digital conversion (ADC) on the error signal peak 
amplitude. The resolution of the ADC operation is A0/4 where A9 is the 
0.3 volt comparator set point. Thus, one data pulse corresponds to a 75 
millivolt increment in the error signal amplitude at the comparator input. 

At the termination of the ADC operation, the logic block sends a reset 
command to the sample-and-compare circuit to return the sample capacitor 
voltage to its quiescent value in preparation for another interrogation 
period. 

During the 6.5 ysec set-up interval, the logic block examines the 
comparator output and then (1) sets the H-switch to route a torque current 
pulse to a dummy load (no torque) or to the torquer, and (2) sets the H- 
switch to route a poritive or negative torque current pulse to the torquer 
if torque is required. At the end of the 6.5 ysec set-up period, the logic 
block transmits a 57. 5 ysec torque command pulse to the H-switch. 

The H-switch torque current source is a precision regulated dc cur- 
rent I (Fig. 5-2). The current I is fixed in amplitude by a 44.353ft pre- 
cision resistor and a precision voltage reference. To change scale factors, 
the voltage reference is switched between 6.25 volt and 1.562 volt levels. 
Therefore, the scale factor can be changed by a factor of four. 

The basic input to the logic block is a 614.4 kHz clock which de- 
termines all the system timing. The other inputs are (1) the two inputs 
from the error signal comparator, and (2) the encode stop signal from the 
sample-and-hold circuit. The outputs are the following: (1) a staircase 
voltage, (2) a 38.4 kHz sync signal to the power supplies, (3) H-switch 
setup signals, (4) torque command signal, (5) gain change command to the 
amplifier and compensator section, (6) scale change command to the bilevel 
dc current regulator, (7) compensation command to the amplifier and com- 
pensator section, (8) sample command to the sample-and-compare circuit, (9) 
a reset command to the sample-and-compare circuit, and (10) the data output 
train. The signals nos. 1, 3, 4, 7, 8, 9, 10 all occur in each 104 ysec 
interrogation interval. The signals nos. 5 and 6 occur only when a scale 
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factor change is demanded. 

The scale decision maker in the logic block commands a scale 
change for the following conditions: the scale factor is switched to 

high if five consecutive torque current pulses are routed to the torquer; 
a switch to low occurs if five consecutive pulses are absent. 

The data train is a sequence of pulses occurring at a 614.4 kHz 
rate. The first pulse occurs at the peak of the 9.6 kHz error signal; it 
is the data sync pulse and is 1 ysec wide. The second pulse contains 
scale factor information — pulse present means high scale, pulse absent 
means low scale. The presence of the third pulse indicates torque applied 
while the absence indicates the converse. The fourth pulse is present for 
positive and absent for negative torque. The next 1 to 31 pulses are data 
pulses, each with 0/4 (75 mv) weighting. Consequently, the error signal 
dynamic range is + 40 or + 1.2 volts. 

ANALYSIS OF SYSTEM BLOCKS 

Logic Blocks 

In this section, the detailed analysis of the operation of the 
Hamilton Standard and the MIT logic blocks will be presented and some com- 
parisons of salient features will be made. 

Hamilton Standard Logic^ . The purpose of the digital function gen- 
erator DFG is to produce in each millisec interrogation period the signals 
shown in Fig. 5-3. The 4.096 MHz clock input signal is counted down to 
produce the 256 kHz pulse train which is used for logic synchronization, 
data generation, mode interrogation signal generation, and basic command 
signal generation. Logic synchronization is accomplished simply by feeding 
the 256 kHz train into the clock inputs of each logic component so that 
logic signal excursions occur on the trailing edge of the sync signal. The 
logic components are standard-line TI SN5400-series devices. 

In the DFG, the 256 kHz pulse train counts down a four-bit binary 
counter (TI SN5493). The 16 kHz output of this counter is used to count 
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Figure 5-3. Signals Generated by the Digital 
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down a second four-bit binary counter. The outputs of the second counter 
feed two binary-to-decimal decoders (TI-SN5442) which provide the 0-16 
blank, 241-256 blank, and a set of gating signals. The gating signals and 
the output signals from the first counter are fed into a set of positive 
nor gates (TI-SN5402) where the gating signals select certain of the first 
counter signals. These selected signals are then transferred through two 

4- bit parallel-in, parallel-out storage registers (TI-SN5495) to the DFG 
output lines. These output signals are the mode interrogation signals PNN, 
PNP, PHN, PHP, PRS, and CRS shown in Fig. 5-3. 

The blanking signals are sent to the quantizer where they limit the 
range of the limit cycle. The time T to the edge of the limit cycle (Fig. 

5- 3) is always at least 1/16 millisec, but never more than 15/16 millisec. 
This assures that the leading and falling edges of the torque pulses are 
never disturbed by limit cycle intrusion. 

The dual scale factor selection (DSFS) logic can be divided into 
four basic sections as shown in Fig. 5-4. They are (1) the scale-factor- 
set J-K flip-flop which is in the "one" state for low scale factor and the 
"zero" state for high scale factor, (2) a nand gate group and J-K flip-flop 
that makes scale factor decisions, (3) a cascade of three four-bit binary 
counters and associated nand gates for providing delay before a scale fac- 
tor change command is transmitted, and (4) a nand gate group for gating 
the 256 kHz data pulses and sync signals into the appropriate output lines. 

Scale factor change requirements are decided in the following manner 
The scale factor set signal (Fig. 5-4), the limit cycle, and the mode in- 
terrogation signals are fed into the "decide scale-factor group". The in- 
terrogation signals PNN and PNP (Fig. 5-3) determine the low mode range; 
the signals PHN and PHP determine the high mode range; the signal PRS sets 
the counter input gates to either count or not count; and the signal CRS 
either counts down the counter or resets the counter to zero. The signal 
CRS is also used as the data sync signal. If the system is in the low mode 
and T < t 2 or T > t^, the counter will be set to count. After 4 counts 
(4 millisec) the scale factor set flip-flop will go to zero and high scale 
factor will be demanded. If the system is in low mode and t 2 < T < t^, 
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the counter will only be reset by signal CRS during each 1 millisec inter- 
rogation period. If the system is in the high mode, and T < t^ or T > t^,, 
the counter will be reset by signal CRS and no scale factor change will be 
demanded. If the system is in the high mode, and t^ < T < t^, the counter 
will be counted down. After 2048 counts (2.048 sec) the scale-factor set 
flip-flop will go to the "one" state and tow scale-factor will be demanded. 
The delay is such that if high mode is needed, the system switches very 
quickly; however, low mode is not set until after considerable delay. 

The scale factor signal (Fig. 5-4) selects the High or Low data 
line and the corresponding sync line in the output of the data gates. The 
limit cycle fixes the length of the 256 kHz data train. The length of the 
data train is always the time T, and an integral number of data pulses is 
always contained in the time T since the edge of the limit cycle is syn- 
chronized with the 256 kHz pulse train that generates the data pulses. 

Note also that the blanking signals in setting the limit cycle range as 
shown in Fig. 5-3 assure that there will always be a minimum of 16 data 
pulses and a maximum of 240 data pulses. 

The results reported above have been experimentally verified on the 
Hamilton Standard logic breadboard constructed at MSFC. The scale factor 
switching was experimentally checked by feeding the logic with a variable- 
width limit-cycle generator which used the 0-16 blanking signal for synchron- 
ization. 

MIT Draper Laboratory Logic* . The clock input to the logic section 
is a 614.4 kHz signal with a 1 ysec wide pulse that occurs every 64th pulse, 
with the other 63 pulses in the overall period having a width of 0.8 usee. 

All inputs to and outputs from the logic section are coupled through dif- 
ferential line receivers and transmitters (NSC DM8820 and DM8830) . These 
differential drivers and receivers retard ground loop problems often en- 
countered in TTL logic, and further permit high level common-mode signals 


*The diagram of this logic block is shown in Fig. 5-2. Circuit details 
are shown in Fig. XI, p. 23, of Reference 2. 
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to exist at the interfaces without disturbing signal processing. Standard 
TTL logic is utilized, with wide application of the newer low-power 74L 
series. Further, considerable use is made of MSI circuits to reduce com- 
ponent count. 

The input clock is fed to a dual monostable RC-controlled multi- 
vibrator, with outputs of 0.8 psec and 1.4 psec wide pulses. The 1.4 psec 
pulse serves as a count-down clock for a set of five J-K flip-flop scalars; 
thus, scaled frequencies of 307, 153, 76, 38, and 19 kHz are available to 
drive a simple (five-bit) staircase generator. Further, the 38 kHz output 
is fed to the drive circuit that controls the DC-to-DC converter power 
supply. The 0.8 psec pulse from the monostable is Nanded with the 614.4 
kHz clock to produce the synchronization signal' (S) and its complement (S). 
These signals are indicated in Fig. 5-5. Either S or S are utilized through 
out the logic section to sync all commands at the 9.6 kHz rate, and to con- 
trol the time at which interrogation of the sample-and-hold circuitry occurs 

If the compensated error signal from the inductive bridge pickoff on 
the gyro is sufficiently low so that the signal to the comparator circuit 
(CA 1 and CA 2 in the SG Processor block) is within the + 300 mv zone, there 
is no torque applied to the gyro and the reference current is routed to the 
dummy load by the R-S logic which turns Q3B "on" and Q3A "off" in the lower 
part of the H-switch. The compensated error signal at the input to the 
comparator circuit is also fed to a JFET sample-and-hold circuit which sam- 
ples the error at the peak. This signal is then compared to the signal 
from the staircase generator, so that with no error signal (zero error sig- 
nal level is 2.50 v) the output from the A to D converter (Data Output Gen- 
erator) is 16 bits. If the error is positive (error signal level at the 
input of the comparators CA1 and CA2 is > +2.50 volts) then the d<:ta output 
is between 17 and 31 bits, while for negative error signal (error output < 
2.50 volts) the data output contains between 1 and 15 bits. Thus, each bit 
is equivalent to 1/4 of the torque pulse quantization for the scale (High 
or Low) in use (which is equivalent to A0/4 * 75 mv, or a data resolution 
of 5.5 arc-sec/4 = 1.37 arc-sec on the low scale). 


DATA 

OUTPUT 


A\ 



Figure 5-5. Control Signals and Data Output in the MIT 
Logic Block. 
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If the signal to the comparator circuit CA1 and CA2 is greater 
than the + 300 mv zone, then output signals A and A (see Fig. 5-2) occur, 
with A in the ”0" state and A in the "1" state for positive error, and 
reversed states for negative error. These signals from CA1 and CA2 drive 
Nand gates which provide the correct torque setup for the Upper and Lower 
parts of the H-switch, and also produce logic signals TM+ and TM- (and 
their complements) to obtain a logic command to order the Compensation 
Signal to be applied to the SG Processor, and further to generate the 
torque signal logic command (TM) to provide the presence of a bit in the 
data train. A real advantage of the ternary system occurs in this sequence, 
since Q^, Q^, Q 2A and Q 2fi are set in their required state before the R-S 
signal is generated that switches the current from the dummy load into the 
gyro torquer winding. Hence, transients occurring during setup of the H- 
switch have no effect. 

When a torque pulse is applied, the logic command TM (and TM) also 
starts a four-bit binary counter (74L93) which acts as the decision maker 
for switching from Low to High torque mode. If five consecutive torque 
pulses are counted, a J-K flip-flop generates the Scale Change Command. 

This command is fed to the data train to establish a "1" indicating High 
Scale, and is coupled to the SG Processor to change the gain by a factor of 
four. An Inhibit pulse is also applied to the R-S switch that inhibits the 
next torque command to allow the H— switch and dc feedback loop to settle. 
Finally, the precision voltage reference (PVR) is switched from 1.6 V to 
6.4 V, or the torquer current is increased by a factor of four. The scal- 
ing logic now continues to monitor the torque signal (TM) and when five 
consecutive torque pulses are absent, the various commands are again given, 
this time to switch back to low quantization. 

As mentioned earlier, the initial four pulses (occupying 6.4 usee) 
in the output data train indicate sync, scale (high or low), torque, and 
sign of the error signal. These pulses are separated from the data pulse 
by an encoding signal which is generated by combining scaled frequencies 
from the Input Scaler along with the 614.4 kHz clock and the output of a 



16-bit binary counter. This timing signal produces a band 5.4 psac wide 
immediately following the sync pulse, within which the encoding informa- 
tion appears in the data train (see Fig. 5-5). 

H-Switch Comparisons 

In the following section, a comparison between the Ternary H-Switch 
of MIT and the Binary H-Switch of Hamilton Standard is made. The H-Switch 
network is one of the most important blocks in the overall electronic re- 
balance loop, since not only the amount of the current through the torquer, 
but also the direction (or polarity) of the current is uniquely controlled 
by the H-switch and its associated error-sensing PVR network. Comparisons 
are made in terms of dc offset currents, plus also transient problems en- 
countered in switching the direction of the current through the torquer 
winding. 

MIT H-Swltch . The schematic of the H-switch section is contained in 

the Howatt*report , and is shown in Fig. 5-6. The torque current ^ tor q ue * s 

precisely controlled by the PVR reference and the negative feedback loop 

composed of R^, A, the Error Amp, and the current-output transistors. Since 

the MIT system is a Ternary design, in the quiescent state (no torque) the 

current I will flow through the dummy load, Q3B, and D7 into R , since 

torque ° P 

the R-S logic will have Q3B "on" and Q3A "off". If torque in one direction 
is required, then logic command signals will turn Q1A and Q2B "on" while 
Q2A and Q1B are turned "off". However, this is an initial setup command 
and occurs before Q3A is turned "on" and Q3B "off". Thus, transient feed- 
through currents in Q1 and Q2 have no appreciable effect. After setup oc- 
curs, logic signals turn Q3A "on" and Q3B "off". The switching transients 
in Q3 will furnish additional energy to R^. However, if the pulse area is 
small, and further if the error-sensing network is slow in relation to the 
transient time, no appreciable transient error should result. In the MIT 
circuit, the dominant open-loop pole in the error-feedback circuit is in 
the input op-amp comparator A (LM 108) and is approximately 1 to 10 Hz. 


•Reference 2, Fig. VI, p. 13, and Sec. 4 (pp. 14-17). 


Driver a 



Figure 5-6. MIT H- Switch 
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Further, diode D7 will help to decouple some of the transient feed-through 
from Rp. Although not Indicated In the schematic for the MIT F-Switch, 
the use of the LM 108A for A does limit the differential input voltage 
(between + and - inputs) to 1 volt, because of the input protective network 
employed in this op-amp. 

In the circuit of Fig. 5-6, not all of the output current I is used 

to drive the torquer, since some of the current feeds the Upper H-Command 

Driver circuit. This latter current splits, with the majority flowing as 

base current drive to the appropriate upper "on" transistor (Q1A or Q1B) 

and the remaining going to circuit ground (The base current does re-enter 

the loop to flow through the torquer). However, it is really not necessary 
• • \ 

that I q be equal to * tor q ue » * n fact, the only real requirement for minimum 

error is that the current through the torquer must all flow through the 
sampling resistor R^. In effect, therefore, a basic theorem for H-switch 
operation is that whatever changes occur abov e the torquer winding creates 
insignificant error in the monitored current loop, whereas changes occurring 
below the torquer winding directly affect system accuracy. 

When the torquer is in the quiescent state (no current flow) , the 
current through R3 and the dummy load is the same as the current through 
the torquer in its active state. 

The offset currents that can exist in the H-switch are those that 
occur in the lower H-switch transistors (Q2A and Q2B) and those that occur 
in the R-S switches (Q3A and Q3B). The model shown in Fig. 5-7 can be uti- 
lized to e inline these offset errors. In this model, it is assumed that 
the requirement is to driv®. the torquer with a current such that Q1A and 
Q2B are "on" while Q1B and Q2A are "off". Also, Q3A must be "on" while Q3B 
must be "off". The resulting direction of current flow through the torquer 
is shown by heavy lines, with other currents indicated by dotted lines. The 
equation for the current sensed by R^ in terms of the actual current through 
the tcrquer 1^ is 

W ■ - twv + + • 


(5-1) 


35. 2 a A (low) 



Flguro 5-7. HIT H- Switch Model for Of f*«t-C\ir)ronc 3»ror An»ly»l». 
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where 1., I., and I. are offset losses, while I and I__ are the base 
14 j dZ dj 

current drive for the "on" state switches (^g and Q^. Since the power 
supplies 5VA and 5VB are floating, these currents should return to the 
supplies and not show up in the output, or Ig^ = Ig2 ' and Igg = Ig^’ . 
Although this is true on a dc basis, it is not necessarily true for 
transient tuvn-on of Q3A. The transient turn-on problem is illustrated in 
Fig. 5-7b, where the lean inductance of the floating power supply re- 
turn has been estimated as 0.1 - 1 yh, and the stray capacitance between 
the floating supply and circuit ground is probably in the 30-100 pf range. 
Since the current Ig^ has two return paths (Ig^ = 1^ + I^.), this means 
that the actual current sensed through R is now I plus I,,, and it 

would thus appear to the feedback network that the current through the 
torquer winding was greater than it actually was. This effect is cer- 
tainly significant under Low Scale conditions, where the torque current 
is 35.2 mA, while the peak value of I „ is approximately 30 mA. Thus, 

D J 

even if only 0.1% of I _ returned to the power supply as the I current , 
this would be a peak transient error of about (0.001) (30mA) / 35 . 2 mA ££ 

850 ppm. It would therefore certainly seem desirable to limit the band- 
width of the error-sensing feedback network for the H-switch to minimize 
transient error. 


Under steady-state conditions, one need only consider the leakage 

current error due to the terms 1^, I and 1^ in relation to * tor q ue * n 

Eq. (5-1). For the 2N2219A transistor used for Q1-Q3 , typical leakages 

at 50°C operating temperature of the torquer would be about 10 nA or less, 

with a temperature coefficient near 1 nA/°C. Thus, the error of these 

terms in comparison to I on the low-scale is 

r torque 


Static Error 


^'torque (low) 


30 x 10 
35.2 x 10" 3 


= 0.9 ppm , 


(5-2) 


and 

Drift Error Ai0.09 ppm/°C 


(5-3) 
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For High-Torque conditions, these numbers would Improve by a factor of 
four. 

Hamilton-Standard H-Switch . A sketch of the Hamilton Standard 
H-Switch is shown in Fig. 5-8, with a more elaborate schematic contained 
in Ref. 5. The essential differences in implementation of the H-switch 
and PVR are (1) that for the Hamilton Standard system the orecision ref- 
erence resistor R is switched (by a factor of 2) whereas in the MIT de- 
P 

sign Rp is fixed while the precision voltage reference is switched, and 
(2) in the MIT H-Switch, the design is based on the use of all bipolar 
transistors with floating power supplies for the lower parts of the H- 
ewitch (as stated earlier, this was necessary to prevent unwanted cur- 
rents from coupling into the sensed current), while in the Hamilton Stan- 
dard circuit, the use of MOSFET pre-drivers for the switching transistors 
obviates the requirement for floating the power supplies. 

A model for the HS circuit is shown in Fig. 5-9, where again the 
main torquer current is shown in heavy lines, with other offset currents 
indicated by dotted liues. Since the drive circuits are isolated from 
the torquer and switches Q5 , Q6, Q9 and Q10 by insulated-gate FETS Q7 , Q8, 
Qll and Q12, there is no quiescent dc gate current present in the torque 
current. There will be transient feedthrough currents through the gate- 
source and gate-drain capacitance of the MOSFETS; however, these should be 
negligible in comparison with the currents encountered in the MIT design. 
It is assumed in Fig. 5-9 that Q5 and Q10 are "on" while Q6 and Q9 are 
"off". The main transient feedthrough will be 1^, I ^ and parts of 1^ and 
1^. The dc offset leakage is due to 1^ and 1^ which gives 


I 

sens. 


I + I, + I- 

torque 4 5 


(5-4) 


From data for the MOSFETS and 2N2219A switches, we obtain for the 
Low Scale torque current of 71.5 mA with maximum offset currents antici- 
pated for both MOSFETS and bipolars of 10 nA, 
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Figure 5 - 8 . Hamilton Standard H- Switch. 








158 


and 


Static Error < 


40 x 10 
71.5 x 10 


0 . 6 ppm , 


(5-5) 


Drift Error £ 0.06 ppm/°C 


(5-6) 


Risetime Effects for the H-Switch . When the torque current is ap- 
plied in the MIT H-Switch, the direction of the current is determined by 
whether the error signal from the gyro pickoff was positive or negative. 

The data shown in the Howatt* report indicates a torque pulse rise-time of 
100 nsec and a fall-time of approximately 300 nsec. Thus, a 1% change in 
rise-time produces a change in torque pulse area of approximately 5 ppm, 
whereas a 1% variation in fall-time produces an area change of about 15 
ppm. These changes are comparable to errors obtained elsewhere in the loop, 
and are therefore of considerable importance. Increasing the rise-and fall- 
times would decrease this error, but may well also increase transient feed- 

j 

through problems. Since some of the same mechanisms control rise-and fall- 
times, there can be partially self-compensating changes produced. 

In the Hamilton Standard H-Switch, since the circuit is operated in 
a Binary Loop, there must be both positive and negative torque pulses ap- 
plied on a continuous basis. Thus, there is an inherent compensation in 
the Hamilton Standard ’^-Switch such that rise- or fall-time variations which 
would decrease the positi e torque area would also decrease the negative pulse 
torque area as well, with a resultant small change in the overall net torque 
pulse area. This compensation will, of course, depend on the matching of 
cross-coupled pairs of drive transistors in the H-switch. 

It would be expected in a low-rate torque environment that imperfect 
compensation in the Hamilton Standard H-Switch would accumulate more torque 
error than would occur in the MIT ternary design. However, in a high-rate 
torque condition, the inherent compensation in the Hamilton Standard design 
should be a significant advantage over the ternary design. 


*Ref . 2, Fig. VIII, p. 17. 
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PV R Circuits 

A comparison is made between the MIT and Hamilton Standard (HS) 
Precision Voltage Reference (PVR) Circuits. Basically, the requirement 
for the PVR is that it produce a scaled precision change in the torque 
current to drive the gyro. The approaches utilized are different for 
these two systems. In the Hamilton Standard PVR a precision resistor 
monitoring network is switched, while the voltage reference is held 
fixed, to produce a change of scale factor. Conversely, for the MIT 
PVR the precision resistor is held fixed, while the voltage reference 
is switched. Numerical values are obtained to enable a comparison be- 
tween drift errors in these two methods. 

Hamilton Standard PVR . A detailed .schematic diagram of the PVR 
is shown in Ref. 3. A simplified model of the PVR is shown in Fig. 5-10. 
The terms suitable for drift error analysis indicated in this figure are 
the following: 

I - the precision current applied to the gyro torquer , 

torque r 

and equal to V REF /R p 

V-p- - the precision Zener diode reference voltage 
REF 

R^ - the precision torque current monitoring resistor 

I - the Zener diode current 

it 

Ip^ - input bias currents for the op-amp 
- the oi :et voltage for the op-amp 

r - the "on" state channel resistance of the FET switches 

on 

r - Zener diode dynamic resistance 

e - the error signal voltage in the feedback circuit 

I - the "off" state channel leakage current of the FET switch, 

off 

The circuit of Fig. 5-10 has the same topological form irregardless of 
whether the system is in a Low or High torque region, hence parameters 
shown in brackets apply to the Low-Scale condition while un-bracketed 


+6.5 V 


3? 



Figure 5-10. Hamilton-Standard Reference, High and Low Scale 
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values Indicate the High-Scale state. Transistor numbers refer to those 
indicated in Ref. 5. 

Analysis of the circuit of Fig. 10 obtains the following equation 
for the error-signal, in the High Scale: 

E - ± V o + V R P + r on« 6) > - h r Z + < V R EF - R p I tcr,ue> 

+ R p [I off (Q1) ' W 9271 + [R p + r on fQ6)] ■ (5 - 7 > 

The drift due to temperature of the error signal e can be estimated by a 
first-order linearization by computing Ae/AT from Eq. (5-7), and using 
values substituted from the data sheets for the appropriate changes. Since 
the ultimate goal is to compare the HS PVR with the MIT circuit, the same 
typical drifts will be utilized for both circuits. Table I below indicates 
the data on drifts obtained from manufacturers data sheets, and values from 
both the HS and MIT reports. 


Table 5-1. 

Drift Data Used for a Comparison Between the Two PVR Circuits 


AV 

“if - ± 5PV/-C 


1^ ■ 75nA(max) 


I 2 " 65nA(max) 


(I^-I^) ~ lOnA(max) 


V REF “ 6 ‘ 35V 


AV, 


~ *+ 5pV/°C 


r ■ 60-ohm 
on 


Ar 

-^ss0.40/*C 

I .. - 2nA 
off 

AI « 

— *0.2nA/°C 
r « 13-ohm 


AR 

- + 3ppm/°C 


R 

P 


44.3530 (MIT) 

310 or 45.50 (HS) 
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Substituting values above and reducing obtains the following 
normalized drift equation for the Hamilton Standard PVR circuit, where 
only terms of significance have been retained. 



(5-8) 


or, substituting values from Table I, 



<8.8 ppm/°C 


(5-9) 


for either High Scale or tow Scale. The long-term drift is also obtained 

similarly using values from the HS report for AV^/At = 23yV/1000 hr, 

AV /At - 10 ppm/ 1000 hr, and AR /At of 12 ppm/1000 hr, gives 
Kbr p 

< 25 ppm/1000 hr . (5-10) 

HIT Draper Lab PVR . The schematic diagram for the PVR is shown in 
the Draper Lab report referred to earlier*. The model utilized for drift 
analysis is shown, for the High-Scale, in Fig. 5-lla and, for the Low-Scale, 
the replacement by the dotted portion of Fig. 5-llb is made. The same para- 
meters are designated as for the HS system, with drift data of Table I uti- 
lized . 

For the High-Scale circuit of Fig. 5-lla, the equation for the error- 
signal is 



*Ref. 2, p. 13, and Appendix I (pp. 32-50). 



Figure 5-11. MIT Reference 
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e - + V + R 

- o p 

- V 

REF 


'We + W* 3 ® ‘ V 

[R 3 + r on (Q1)][I l ' W" 2 ” > 


(5-11) 


where R^(7.5K) and R2(2.5K) are precision resistors used to form a 4:1 vol- 
tage divider, and R^ is a resistor chosen to offer the same impedance to the 
op-amp, regardless of the scale employed. 

The computation of Ae/AT from Eq. (5-11) is obtained as was that for 
the HS circuit. Substituting the values from Table I, and reducing, the 
following equation for drift is obtained for the MIT PVR circuit: 



Ae 

AT 

High 


+ 



AR 

( — E) 
'‘AT ' 


AI AI (Q2) 

[ r 3 + r on (Q 1 )H'^f - + ] 


(5-12) 


Substituting values from Table I, 


9ppmrc 


(5-13) 


The long-term drift is similarly, 


-U < 25 ppm/1000 hr + ~ [R^ + r (Ql)] x 

V™ ) Atime "* V_ Dt , 3 on 


REF, 


REF 


All AI off (Q2) 

[ Ar~ + "At J * 


(5-14) 
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For Low-Scale, using the circuit of Fig. lib produces (with k 


R 2 /(R x + R 2 )J 


V REF 


, [ av 

1 J o | 

~ + c + 


Ac 

AT “ - k V 


k AV 


REF 


REF 


AT AT 


k V REF AR n 

+ it ( AT^ + [k R 1 + r on (Q2)1 x 


m, m oK « 13) ak 

f — - 2±± 1 4- V (— ) 

1 &T AT 1 REF '‘AT' 


(5-15) 


or 


Temp, drift (Low-Scale) £,14.4 ppm/°C . (5-16) 


A similar equation for long-term drift for the Low-Scale can be obtained 
and is similar to Eq. (5-15). 

The basic conclusion to be obtained from this analysis is that al- 
though the drift races are comparable for High-Scale torque currents, under 
Low-Scale conditions it is better to leave the reference comparison voltage 
fixed and switch the 'precision resistor R^. It is also apparent that switch- 
ing the reference comparison voltage allrws other drift terms in the circuit 
to be important, as evidenced by comparing the terms of Eqs. (5-8) and (5-15). 

No comparison has been made of the switching transients involved in 
changing scale factor. However, it is apparent that in switching either re- 
sistors or voltage references, transient feed-through errors will result, 
unless the loop is given sufficient time to settle before torque is applied 
to the gyro. This can be done quite easily for the Ternary system, but not 
for the Binary loop. Further, the op-amp (A in Figs. 5-10 and 5-11) must 
have sufficient differential input rating so the switching transients do not 
cause junction breakdown. 
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Noise Considerations 

The effect of electronic noise on system resolution will be con- 
sidered in this section. The calculation details for obtaining a noise 
figure-of-merit for both the Hamilton Standard and MIT Systems will be 
presented and the noise performance of the two systems will be compared. 

MIT System . The electronic noise contribution to the system un- 
certainty can be determined by computing the rms noise voltage at the in- 
put to the comparator block of Fig. 5-2. Since a voltage increment at the 
comparator input represents a known angular increment in gyro float position, 
the noise can be expressed as rms arc-seconds of uncertainty. At the com- 
parator input, a 300 mv. voltage increment corresponds to 5.5 arc-seconds 
on the low scale and 22 arc-seconds on the high scale. This A9 of 300 mv. 
is the weight of one torque pulse. 

The white noise voltage at the comparator input is given by 

e = A_. |4kTR Af "I 1/2 , (5-17) 

where k is Boltzmann's constant, is the voltage gain from the amplifier 
input to the comparator input, T is the absolute temperature, R^ is the 
equivalent input noise resistance of the amplifier and compensator block 
(Fig. 5-2), and Af is the noise bandwidth from amplifier input to compara- 
tor input. Since the gain of the input stage of the amplifier is about 20, 
the noise should be controlled by this input stage. Assuming the error- 
detector output impedance is near zero, the equivalent input noise resis- 
tance should be about 8 kilohms (7 kilohms attributed to the input network 
and 1 kilohm attributed to the UC 4250 input op-amp) . The dominant low- 
frequency 3 db point is about 3.4 kHz and is controlled by the coupling 
network between the second and third stages of the amplifier. The domi- 
nant high-frequency 3 db point is set by the input network and is about 8 
kHz. Thus, a reasonable approximation to the noise bandwidth is 5 kHz. 

The gain is about 4000 on the low scale. Using these values in Eq. 

(5-17) gives 
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e =3.3 millivolts rms. (5-18) 

nw 

Therefore, the noise referred to the gyro float motion is 

„ w 5.5 arc-sec s 

A0 = (3.3 mv rms) ( — ) , 

n 300 mv ’ 

or 

A0 =0.06 arc-sec rms. (5-19) 

n 

This value of A0^ should be considered a minimum value because low-frequency 
electronic noise and electronic noise generated by the error signal detector 
have been neglected; however, these two additional sources would not be ex- 
pected to increase the noise of Eq. (5-19) by more than about 20%. 

Hamilton Standard System . The sensitivity of the error-signal detec- 
tor in the Ham. Std. System is given as 1.39 millivolts per arc-sec*. The 
white noise f igure-of-merit can therefore be calculated by first finding 
the equivalent noise ^oltage at the preamplifier input (Fig. 5-1). 

The IC^ stage of the preamplifier will make a significant contribu- 
tion to the input noise, since it is preceded by a unity-gain section. The 
IC^ gain-of-6 should assure that only the preamplifier need be considered 
in the noise calculation. 

The estimated total input-noise resistance R^ is composed of the 
following contributions: (1) differential source follower, 2K; (2) differ- 
ential emitter follower, 2K; (3) differential inverter, 4K; and IC^ circuit, 
20K. These noise components sum to a value of 28K and are estimated to be 
caused primarily by circuit resistors. 

Before a final value of R is determined, the noise contribution of 

n 

the 1-meg gate resistors in the input source follower must be computed. The 
output ring capacitance of the error signal detector (a capacitance bridge) 
is 116 pf.* The 3db frequency associated with the gate resistance is 


*Ham. Std. Report, Ref. 3, p. 125. 
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0 159 

f = rpr * 685 Hz . (5-20) 

8 (2 x 10 ) (1.16 x 10 ) 


the effective noise resistance at 50 kHz then becomes 


R 

ng 



< 6 - 85 * <2 - 10 6 > 

(5 x 1CV 


(5-21) 


or 

R *= 360 ohms 
ng 

Thus, the contribution of the gate resistors is negligible. 

Taking a value of R^ - 28K for the total white noise resistance, the 
input equivalent-noise-voltage density is 

ENV = AkT R = 21.7 nanovolts/Hz^^ . (5-22) 

n 

The noise bandwidth of the 50 kHz filter is about 5 kHz, thus, the equi- 
valent noise voltage is approximately 

ENV t = 21.7 /5 , 000 * 1.54 microvolts rms . (5-23) 

Referring this to the detector signal gives 

A0 ** (1.54 x 10 -6 ) r - 0.0011 arc-sec rms . (5-24) 

n 1.39 x 10' J 

Comparisons . The electronic noise of the Ham. Std. rebalance elec- 
tronics is much less then that of the MIT system. Yet, the equivalent noise 
resistance of the MIT system is less than that of the Ham. Std. system. The 
reason for the poorer signal-to-noise ratio of the MIT system is the much 
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lower sensitivity of the error-signal detector. This lower sensitivity is 
due not only to basic differences in the detectors (MIT is an inductance 
sensor; Ham. Std. is a capacitance sensor) but also to the higher frequency 
(50 kHz vs. 9.6 kHz) error signal employed by Ham. Std. Undoubtedly, the 
MIT detector could be made much more sensitive with a consequent system 
signal- to-noise ratio comparable to that of the Ham. Std. System. 

The higher frequency error signal gives the Ham. Std. System a sig- 
nificant advantage if low-frequency electronic noise is a problem. For 
practically all modern semiconductor devices, a frequency of 50 kHz is well 
above the frequency range where low-frequency noise is important. Semi- 
conductor devices must be selected to have negligible low-frequency noise at 
9.6 kHz. 

Neither the Ham. Std. nor the MIT preamplifier is a low-noise design. 
However, improving the noise performance of either preamplifier without 
sacrificing common-mode rejection and impedance requirements will not be a 
trivial exercise. 

Common-Mode Rejection 

Excellent common-mode rejection in the preamplifiers of both systems 
is essential for eliminating common-mode signals that could contaminate the 
error signal. Common-mode signals of several volts at the preamplifier input 
should have less effect on the preamplifier output than the millivolt-level 
differential signals from the detectors. Consequently, a common-mode re- 
jection CMR of no less than 80 db would be indicated for the preamplifiers 
used in both the MIT and the Ham. Std. Systems. 

No extensive detailed analysis has been completed on either system's 
preamplifier; however, preliminary results indicate that the CMR is probably 
adequate. 

The Ham. Std. System has an elaborately designed preamplifier using 
discrete transistors and a monolithic differential amplifier with great em- 
phasis placed on high CMR. The achievement of high CMR in the Ham. Std. 
System is greatly enhanced by having low-gain (about 6) before the bandpass 
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filter and by using a 50 kHz error signa] which should be much higher in 
frequency than most expected common-mode signals. 

The MIT preamplifier is a single monolithic differential amplifier 
which could possibly be selected to have a CMR of about 80 db at 9.6 kHz. 

The lower frequency bandpass (about 3 kHz to 8 kHz) of the MIT amplifier 
probably opens up this system more to CM signals; however, the much lower 
senjitivity of the error signal detector and the fact that monolithic op- 
erational amplifiers have maximum common-mode rejection at low frequencies 
may be partially compensating. 

Compensation 

The methods of compensating the rebalance loop used in the Ham. Std. 
and the MIT Systems are very different. In this section, the two types of 
compensation are described and their salient differences delineated. 

Hamilton Standard Rebalance Loop Compensation . In the Ham. Std. Sys- 
tem, a cascade compensation network was used to (1) provide a stable loop 
with sufficient phase and gain margins, (2) discriminate against gyro float 
angular vibrations that tend to saturate the loop amplifiers, and (3) at- 
tenuate signals arising from the gyro wheel rotation. For the gyro used in 

3 

th_ Ham. Std report the transfer function in the frequency domain for the 
gyro alone was approximately 

K 

C(j f ) = 8 f , (5-25) 

Jf(1 + J 1050 ) 

where f is the frequency in Hz. The transfer function for the compensation 
network (Fig. 5-1) used with this gyro was 

K (1 + jy^yy) 

C (jf ) . (5-26) 

(1 + ^0.016^ ^ + JlOO 5 

Except for some poles at high frequencies that should be negligible, 
the overall frequency domain transfer function viewed by the torque current 
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pulses is the product -C(jf) G(jf). Clearly, the overall transfer function 
can not be approximated by that of a pure integrator. It appears that making 
the loop approximate a pure integrator was not a design goal for the Ham. 

Std. System as it was for the MIT System. 


It may be difficult to achieve system stability for very high loop 
gains since the maximum phase shift of the -C(jf) G(jf) product is 450°. At 
a frequency of 5 kHz, where the bandpass filter begins to contribute phase 
shift, the frequency dependent part of C(jf) will have attenuated the loop 
transmission by a factor of about 3500. This attenuation in addition to 
that of the gyro dynamics should assure that the phase shifts of the high- 
frequency poles in the rebalance electronics will not cause system instability. 

Even though the Ham. Std. compensator transfer function is more complex 
than that of the MIT compensator, the actual implementation can be rather 
easily achieved using standard operational amplifier networks. 


MIT Rebalance Loop Compensation . The MIT compensator ic simple and 
is designed to cause the effective transfer function from the torquer input 
to the comparator input (see Fig. 5-2) to approach that of a pure integrator. 
The detailed development of this method of compensation is presented in the 
M.S. Thesis by Lory^. The compensation technique is shown in Fig. 5-12. 


The overall transfer function from point T to Point C in Fig. 5-12 
is given by 


T(s) 


K + H 2 (s) [sCt^s + 1)(t 2 s + 1)] 
s(t 1 s + 1) (x 2 s + 1) 


(5-27) 


Usually the gyro time constant x^ is much greater than the torquer 
time constant x 2 (x^ is typically greater than 100 microseconds). For this 
usual case T(s) can be made equal to 1/s if (from Eq. (5-27)) 


H 2 (s) 


1 

x^s + 1 


+ 


1 - K 
s(x^s + 1) 


(5-28) 


The first term on the right side of Eq. (5-28) is a simple RC low-pass 
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network. The second term is a simple RC low-pass network cascaded with an 
integrator. The two terras can be realized, therefore, as the parallel com- 
bination of simple networks. 

It is especially significant to note that the compensation only need 
be applied during the actual torquing and can be driven by any rectangular 
pulse with edges coincident with the actual torque current pulse and ampli- 
tude directly proportional to the torque command pulse generated in the logic. 
(Torque current pulse changes are accompanied by compensating loop-gain 
changes to hold the overall loop-gain constant.) 

In the actual implementation of the MIT System, the compensation com- 
mand generator is driven by the torque command pulse. The compensation net- 
work is a simple RC low pass circuit which provides accurate compensation 
only for the case K ■ 1. 

It may appear a bit disturbing that the MIT compensation technique is 
equivalent to actually feeding a signal that is proportional to the torque 
current pulse around the gyro, detector, and error-signal amplifier. How- 
ever, since the overall transfer function is nearly a pure integrator, the 
signal that appears at the comparator input is directly proportional to the 
area of the torque current pulses. Now, since the float angle increment A0 
is directly proportional to the area of the torque current: pulses, then the 
signal at the comparator input is directly proportional to A0 . Thus, the 
sample-and-compare circuit will indeed generate data pulses proportional to 
A6 which, of course, is the desired result. 

The ease of implementation, the conceptual simplicity and the expected 
reliable performance of the MIT compensation technique mak i the technique an 
attractive way of compensating ternary systems with the consequent elimination 
of pulse bursting. One possible disadvantage of this compensation technique 
is the dependency of ^(s) on the gain K. If after ^(s) is implemented, 
either the gain of ^(s) or the gain of the preamplifier changes, then the 
compensation is no longer perfect. We have not yet determined the detailed 
effects of such gain drifts. 
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Another important consideration in compensation design is the effect 
of the width of the torque pulse on the system response. Since the MIT 
torque current pulse is only 97.5 microseconds compared to a maximum width 
of 15/16 milliseconds for the Ham. Std. System, it is much easier to make 
the loop appear to be an integrator for the MIT torque pulse. When the 
torque current pulse is sufficiently short relative to the gyro dominant 
time-constant (the pseudo-integrator), then several error-signal processor 
design constraints are relaxed. Among the relaxed constraints are those 
affecting the dynamic range and the spurious signal discrimination circuits 
(for instance, gyro float oscillations at the forced limit-cycle frequency) 
in the error-signal processor. 

OVERALL SYSTEM CONSIDERATIONS 

Three vital considerations in the comparison of all rebalance elec- 
tronic schemes are resolution, accuracy, and reliability. Resolution and 
accuracy are usually interdependent and must be considered simultaneously 
in the design of rebalance electronics. For example, improvement in reso- 
lution at the expense of accuracy could result in a degradation of overall 
system performance. Also, reliability is obviously closely connected since 
reliability involves many alternatives not only in technique but also in 
selection of components. These choices of techniques and components must 
always be made carefully. 

In this section, the resolution, accuracy, and reliability of the 
Ham. Std. and MIT rebalance electronics will be discussed and comparisons 
made where possible. The presentation will draw heavily on the results of 
the analyses given in the previous sections. 

Resolution 

To compare torque current-pulse area quantization and data pulse 
quantization, the gain of the system comprising the torquer , gyro, and 
error-signal detector can be removed by considering torque current-pulse 
area increments only at the input to the torquer. Inis normalization al- 
lows a direct comparison of the inherent quantization produced by the 
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rebalance electronics only. 


The Ham. Std. System has a one-millisecond torque current pulse with 

both positive and negative area. The pulse can be changed in 3.9 psec 
1 -3 

(ttt x 10 sec.) increments. Increasing the positive current width by 3.9 
zoo 

psec also decreases the negative current width by 3.9 psec. Thus, the 
smallest area increment AA is given by 


AA - 3.91 - 3.9 (-1) 


“h- 


7.81 nanoarcp- seconds 


(5-29) 


where I is the amplitude of the torque current pulse in milliamperes . This 
increment in current-pulse area will produce a one pulse increment in the 
data train; thus, the data quantization AD^ is also 

AD = 7.81 nanoamp-seconds . (5-30) 

H 


The finest quantization occurs in low-scale operation where I is 71.5 ma. 
Hence, the data quantization for the Ham. Std. System is 


AD = (7.8) (71.5) = 558 nanoarap-seconds , 

n 


and also 


“h 


= 558 nanoamp-seconds 


The HIT torque current-pulse width is constant and is 97.5 micro- 
seconds. The torque current pulse quantization is therefore 


AAjj - 97.51 nanoamp-seconds , (5-31) 

If I is the amplitude of the current pulse in milliamperes. The data pulse 
weighting is such that four data pulses correspond to one torque pulse. 
Hence 

AD„ = 24.41 
M 


(5-32) 
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The low-scale current amplitude is 35.2 ma; consequently, the numerical 
value of the MIT System data quantization is 

AD__ = (24.4) (35.2) = 860 nanoamp-seconds 

M 

The difference in data quantization for the two systems is not very 

great - AD = 1.54 AD - but it is also not insignificant. There is an 
M H 

appreciable difference in torque current-pulse quantization - AA^ = 6.2 AA^. 

In conclusion, then, the data resolution of the Hamilton Standard System is 
somewhat higher than that of the MIT System, while the torque current-pulse 
resolution of the Ham. Std. System is much higher. 

Another important resolution consideration for some strapdown gyro 
applications is how precisely the position of the gyro float angle is known 
as a function of time (time resolution). The MIT System has an apparent ad- 
vantage in this case since the error signal is sampled at a well defined time 
and at a rate that is about an order of magnitude higher than in the Ham. 

Std. System. There is, however, some cause for concern if the MIT compensation 
is not perfect; that is, either too much (overcompensation) or too little 
(undercompensation) signal is fed through the parallel compensation feedback 
path into the comparator input. This imperfect compensation, which could arise 
from component change due to aging or temperature variations, may possibly 
cause short term error in data output with a consequent degradation in the 
time resolution of the system. If the compensator parameters drifted very 
much, perhaps new pulse bursting modes may be introduced. Further study of 
the MIT compensation technique is needed to provide definitive conclusions 
concerning the effects of imperfect compensation on time resolution of the 
system. 

Accuracy 

Assuming that the torquer has negligible non-linearity, the accuracy 
of both the MIT and the Ham. Std. Systems should depend primarily on the 
stability of the torque current-pulse amplitude regulator and on the perfor- 
mance of the H-switch. There are two secondary sources of inaccuracy that 
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could also be of significance. These are (1) instabilities in the analog- 
to-digital conversion for data generation and (2) spurious signal contami- 
nation of the error signal. 

We will discuss the four sources of inaccuracy identified above, 
starting first with the secondary sources. 

Analog-to-Digital Conversion . The analog-to-digital conversion, 

ADC, for data generation in the Ham. Std. System is accomplished by detect- 
ing the time t that a ramp signal crosses a slowly varying voltage level 
proportional to tie error signal. The time t is sensitive to gain drifts 
in the error signal amplifiers and to the slope of the ramp. However, the 
entire ADC operation is accomplished inside the, system loop so that for 
practical values of loop transmission, these inaccuracies should be made 
negligible by the feedback mechanism. 

A very different i.ituation exists in the MIT System. The ADC is 
connected parallel to the loop and is, therefore, not included in the loop. 
Thus, the stability of the calibration of the data output depends on the 
stability of the staircase generator. The critical resistors in the stair- 
case generator should have stability comparable to that of the current sen- 
sing resistors in the torque current regulator. We are presently examining 
ways of eliminating this deficiency in the ternary system. 

Spurious Signals . An important spurious signal is electronic noise 
superimposed on signals at points in the system where torquing decisions 
are made. Since this noise will be gaussian distributed about a zero mean 
value, its effect on system accuracy should be averaged to very small values 
over several torque pulses. Because of this averaging, the rms noise errors 
computed in the previous section on Noise Considerations are pessimistic. 

The inaccuracies due to other spurious signals that are not corre- 
lated with the error signal should also be averaged over several torque 
pulses to negligible values. Correlated signals may possibly produce ser- 
ious inaccuracies; thus, care in layout and attention to ground-loop elimi- 
nation is essential in implementation of rebalance systems. 
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Effects of spurious signals on accuracy must usually be determined 
experimentally . 

Current Regulator . The current regulator circuit is composed of the 
dc feedback loop containing the H- Switch, the Precision Voltage Reference, 
the comparator op-amp A and the error amplifier. This circuit was shown in 
Fig. 5-2, and in more detail in Figs. 5-6 and 5-8 for the two systems. It 
is absolutely essential for system accuracy that the weight of each torque 
pulse be precisely known; i.e., that the amplitude of the current pulse and 
its pulse width be accurately maintained, or more correctly that the area of 
the pulse be constant. The pulse width is basically controlled by the clock 
accuracy, and should contribute the same inaccuracies whether the system were 
binary or ternary. Calculations made earlier allow some comparison of pulse 
rise- and f all-time inaccuracy, and also inaccuracies in current amplitude. 

A worst-case comparison of the MIT and Ham. Std. current regulator designs 
can be made by adding these various inaccuracies, or 

Overall Accuracy (worst-case analysis) 

& error (PVR) + error (Sensed Current) 

+ error (rise- and fall-times) (5-33) 

Collecting results gives for the Low-Scale mode the following comparisons: 

Acc. (worst-case) < (14.4ppm/°C and > 25ppm/1000 hr) 

Mil — 

+ (0.9ppm and 0.09ppm/°C, plus transient effects) 
+ (JKlOppm for a 1% change), (5-34) 


and 

Acc. (worst-case) < (8.8ppm/°C and 25ppm/1000 hr) 

nb — 

+ (0.6ppm and 0.06ppm/°C) 


+ (<10ppm for a 1% change) 


(5-35) 
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Reliability 

No comparison has been made between the MIT and the Ham. Std. Systems 
regarding reliability, as this would involve considerable in-depth experi- 
mental analysis of the two. Qualitatively, there seems to be no reason one 
circuit should be any more reliable than the other. The considerable use 
of MOS transistors in the Ham. Std. design would cause problems if the cir- 
cuits operated for an extended time in a radiation zone, since it has been 
found that the threshold voltages of most enhancement-mode MOSFETs in- 
crease considerable in a radiation field. 

SUMMARY AND RECOMMENDATIONS 

The results of our analyses of the rebalance electronics of the two 
strapdown gyroscopes (MIT ternary and Ham. Std. binary) have not led to a 
conclusion that one system is clearly superior to the other. We have used 
our analytical results to produce a tabular comparison (Table II) of the 
two systems in seventeen important areas. Each system has a significant 
set of both relative advantages and disadvantages. We cannot distinguish 
an appreciable difference in the overall quality of electronic engineering 
design principles applied in the development of the two systems. 

We recommend that the feasibility be explored of developing a system 
that would incorporate the best features of both the ternary and binary 
systems studied. Such a system perhaps would have the capability of opera- 
tion in both ternary and width-modulated binary modes with built-in decision- 
making apparatus for causing the system to operate in the mode that is best 
for the rate environment in which the gyroscope finds itself. For example, 
a zero rate condition would probably require a ternary mode to eliminate 
cumulative torque-current bias errors; yet, a high rate environment may 
need a binary mode to take advantage of the inherent rise- and fall-time 
compensation in the binary H-switch. 

Early in the development of a new system, an in-depth experimental 
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study of the two types of H-switch (binary and ternary) must be made since 
major design decisions will depend strongly on fundamental limitations in 
the performance of the H-switches. Also, of central importance are the 
loop transmission (which strongly influences system accuracy) system sta- 
bility trade-offs and the design compromises among high interrogation rates, 
electronic implementation limitations, time resolution of data output, and 
accuracy. 
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CHAPTER 6 

PROBLEMS FOR FURTHER STUDY 

Continued effort is needed to further develop the technology for 
the high reliability strapdown systems using TDF gyros. In particular, 
the following activities are recommended: 

1. Development of optimum redundancy configurations for 
systems using more than six gyros. 

2. Analysis of the system performance of the rebalance 
loops of a TDF sensor. 

3. Further development of the sensor performance management 
concept to see if any more sophisticated detection and 
identification schemes can be employed to improve the 
overall system performance. 

In the area of rebalance electronics, more work needs to be done 
to arrive at a design which is acceptable for NASA’s applications. The 
following activities are desired. 

1. Completion of the evaluation of MIT's and Hamilton 
Standard's electronic designs. 

2. With the help of the background gained by evaluating 
the above designs, develop a better design of the re- 
balance electronic circuit. 

3. Review of all applicable technology in the area of 
pulse rebalance electronics. 

At present, various concepts are available for using SDF gyros as 
strapdown sensors. A study and evaluation effort is needed to extract and 
modify the best concepts so they can be incorporated together ,_o form the 
best overall system. The following questions need to be answered. 

1. Which gives a better system performance between the 
binary and the ternary pulsing concepts? 
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2. Among the known concepts of failure detection and 
faulty sensor identification, which one is the best? 

3. What are the optimum sensor arrangements when more 
than six SDF sensors are used? 

Can any better concepts be developed? 


4. 



